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1.  FOREWARD 


This  is  the  final  rep)ort  arising  from  our  research  program  involving  the 
investigation  of  processing/microstructure  relationships  in  selected  ceramic  systems.  The  first 
year  of  this  program  was  dedicated  to  the  development  and  detailed  assessment  of  the  novel 
neutron  scattering  and  x-ray  diffraction  techniques  required  for  the  measurement  of  ceramic 
microstructure  evolution  as  a  function  of  sintering  and  chemistry.  The  next  two  years  of  the 
program  were  used  to  study  several  systems  in  detail,  and  has  led  to  some  remarkable 
conclusions  which  have  changed  the  way  in  which  the  sintering  stages  in  different  systems  are 
understood. 

The  sintering  behavior  of  a  ceramic  body  and  the  properties  of  the  ceramic 
product  depend  directly  on  the  internal  microstructure.  The  goal  of  our  research  program  with 
the  ARO  has  been  to  investigate  well-characterized  systems  and  to  measure  microstructure 
evolution  during  ceramic  processing  in  order  to  obtain  an  improved  understanding  of  the 
relationships  between  processing  and  microstructure.  This  approach  has  led  to  progress  in 
improving  process  models  and  to  improved  predictability  of  product  microstructure. 
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2. 


INTRODUCTION 


Starting  with  the  multiple  small-angle  neutron  scattering  formalism  originally 
developed  by  Berk  and  Hardman-Rhyne*  at  NIST,  we  applied  it  to  such  questions  as: 

o  What  effect  does  the  green-body  density  have  on  pore-size 
evolution  and  grain-size  evolution  during  sintering? 

o  What  is  the  role  of  a  sintering  aid  in  the  intermediate  and 
final  stages  of  sintering? 

o  What  signals  the  beginning  of  final  stage  sintering  and  how 
does  the  porosity  evolve  as  full  theoretical  density  is 
approached? 

o  How  do  the  sintering  mechanisms  (surfece  and  volume 
diffusion;  viscous  flow)  change  the  porous-phase  evolution? 

The  answer  to  this  last  question,  for  example,  is  beginning  to  have  a  major  impact  on  the 
community  because  the  results  differ  significantly  from  classical  well-accepted  theory.  The 
answers  to  all  of  these  questions  will  be  described  in  the  neutron  scattering  section  (Section  3) 
which  follows. 

High  temperature  x-ray  diffraction,  incorporating  several  ceramic-specific 
features,  has  been  used  to  investigate  in-situ  sintering  in  terms  of  the  evolving  grain  sizes  and 
phases  of  several  systems.  In  addition,  we  have  studied  the  sintering  kinetics  of  materials 
synthesized  at  NIST  under  the  support  of  this  program.  Finally,  the  evolution  of  pore  sizes 
and  size  distributions  during  final-stage  sintering  has  been  investigated  using  small-angle  x- 
ray  scattering.  The  latter  results  indicate  that  for  many  applications,  a  product  with  less  than 
full  theoretical  density  leads  to  properties  superior  to  a  product  with  nearly  full  density.  These 
results  will  be  described  in  the  x-ray  section  (Section  4)  which  follows. 
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3. 


SMALL-ANGLE  NEUTRON  SCATTERING  CHARACTERIZATION 


The  central  thrust  of  this  research  program  was  the  application  of  advanced 
scattering  techniques  to  study  the  processing/microstructure  relationships  in  real  ceramic 
materials,  leading  to  an  improved  understandi»ig  of  ceramic  process  models.  In  this  study,  it 
was  possible  for  the  first  time  to  follow  particle  sizes,  pore  sizes,  size  distributions  and 
morphology  from  the  earliest  "green"  state,  through  intermediate  stage  sintering,  through  to  the 
final  stages  of  sintering. 

The  first  system  that  we  studied  was  undoped  o'-alumina,  which  served  as  a 
model  material.  Our  intention  was  to  test  the  quantitative  capabilities  of  our  new  scattering 
techniques  with  a  system  which  is  well  understood.  We  measured  the  effective  pore  radii  for 
alumina  56-98%  dense  and  found  that  we  could  not  only  quantify  the  pore  sizes  but  we  could 
also  identify  the  transition  from  intermediate-stage  sintering  to  final-stage  sintering  from  the 
morphological  change  which  takes  place  in  the  shape  of  the  pores  present  at  each  stage.  In  the 
intermediate  stage  we  measure  the  cross  sections  of  the  interconnected  (i.e.,  open)  pore 
channels,  whereas  in  the  final  stage  we  measure  the  diameters  of  isolated  pores.  Our  ability 
to  measure  the  density  at  which  the  transition  from  intermediate  to  final  stage  takes  place 
turned  out  later  in  this  program  to  be  essential  for  imderstanding  the  effect  of  sintering  aids 
on  microstructure  and  the  effect  of  green-body  density  on  product  microstructure.  These 
results  will  be  detailed  below.  At  this  time,  we  show  in  Figure  1  the  effective  pore  radius  as 
a  function  of  sintering  for  undoped  alumina.  For  densities  between  55%  and  87%  of 
theoretical  density,  the  effective  pore  radius  remains  nearly  constant.  This,  at  first  surprising, 
result  for  intermediate  stage  sintering  is  actually  consistent  with  sintering  modules  in  which  the 
stable  topological  decay  of  an  interconnected  pore  network  is  taking  place.  This  constant- 
radius  result  is  confirmed  by  mercury  porosimetry  results  and  Porod  scattering  results  on  the 
same  samples.  For  densities  between  90%  and  98%  dense,  the  radii  of  the  remaining  pores 
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increases  monotonically.  We  have  identified  these  radii  as  being  a  measure  of  the  pore  size 
distributions  of  the  final  stage  of  sintering,  in  which  the  pores  are  isolated  from  one  another 
and  the  decay  in  percent  porosity  is  no  longer  stable.  Thus  MSANS  and  Porod  scattering  have 
identified  for  the  first  time  where  the  transition  to  isolated  porosity  takes  place,  and  has  led  to 
strong  confirmation  of  the  topological  decay  model^  for  intermediate-stage  sintering,  which  until 
now  has  probably  not  received  the  attention  it  deserves. 


Fig.  1 .  Effective  pore  radii  for  the  13>series  alumina  as  a  function 
of %TD  derived  from  MSANS,  from  Porod  scattering,  and  from  MIP. 

Next  we  studied  the  effect  of  green  density  and  the  role  of  a  sintering  aid  on  the 
densification  during  sintering  of  alumina.  It  is  well  known  that  high  green-state  porosity  leads 
to  difficulties  in  sintering,  and  it  is  also  well  known  the  addition  of  small  amount  (0.25  wt  %) 
of  MgO  leads  to  a  much  finer-grained  and  denser  ceramic  product.  It  was  not  known, 
however,  why  a  low  green-body  density  could  not  be  overcome  by  sintering  or  the  mechanism 
by  which  the  addition  of  MgO  inhibits  grain  growth  to  form  a  much  finer-grained  denser 
ceramic  product  with  fewer  incipient  porous  flaws.  In  research  supported  by  this  program, 
we  examined  the  sintering  of  alumina  with  and  without  additive,  and  we  followed  the  course 
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of  the  pore  size  distributions  starting  with  green  densities  from  46%  to  53%.  The  results 
showed  unequivocally  that  the  initial  connectivity  in  the  green  body  detcTOined  the  size  of  the 
pores  in  the  intermediate  sintering  stage,  which  are  always  constant  (see  Figure  2).  Further, 
it  was  demonstrated  that  the  density  at  which  the  transition  from  stable  intermediate-stage 
sintering  to  final-stage  sintering  takes  place  is  determined  in  part  by  the  green  body  density. 
The  denser  the  green  body,  the  later  the  transition  takes  place,  and  the  smaller  the  isolated 
pores  which  remain  at  the  end  of  the  final  stage. 


2.  Effective  pore  radii,  as  a  function  of  %TD  for  all 

four  series  of  samples.  The  number  appearing  near  the  leftmost 
symbols  is  the  average  green  body  density  for  each  series. 

With  regard  to  the  effect  of  sintering  aid,  it  was  learned  that  the  addition  of 
0.25  wt  %  MgO  can  overcome  some  of  the  problems  associated  with  low  green  body  density. 
Our  MSANS  data  indicated  that  even  with  low  green  body  density,  the  onset  of  final  stage 
sintering  is  delayed.  This  leads  to  finer  grains  and  smaller  pores  in  the  product  material 
because  the  stable  intermediate  densification  stage  is  extended  to  densities  well  over  92%  of 
theoretical. 
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Although  it  is  difficult  to  achieve  a  concensus  regarding  many  aspects  of 
sintering  mechanisms,  there  is  general  agreement  that  the  driving  force  for  sintering  is  the 
reduction  in  surface  free  energy.  The  shape  of  the  minimum  free  energy  surface  influences 
strongly  the  microstructure  evolutions  during  processing.  In  another  research  problem 
investigated  with  the  support  of  this  ARO  program,  we  measured  the  pore  sizes,  size 
distributions,  surface  areas  and  number  distributions  of  the  porous  phase  in  glassy  silica  and  in 
polycrystalline  alumina.  In  this  way  we  had  a  probe  of  two  very  different  sintering 
mechanisms:  viscous  flow  in  the  processing  of  the  glass,  and  surface  and  volume  diffusion  in 
the  processing  of  the  polycrystalline  material.  The  goal  was  to  gain  a  quantitative  measure  of 
the  structural  evolution  when  different  sintering  mechanisms  dominate  densification.  The 
remarkable  result  of  this  study  was  that,  from  a  topological  point  of  view,  both  systems  evolve 
the  same  way.  This  is  shown  in  Figure  3,  where  we  show  the  pore  volume  fraction,  0, 


50  60  70  60  90  100 


%  TD 

Fig.  3.  Pore  volume  fraction,  normalized  to  S/V,  as  a  function  of 
densification  (%  TD)  for  silica  and  alumina. 
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normalized  to  the  pore  surface  area  per  unit  volume,  S/V,  as  a  function  of  densification.  The 
points  for  both  systems  for  this  topological  parameter,  0/(S/V),  fall  on  a  universal  curve. 
However,  in  terms  of  the  actual  pore  sizes  and  morphology,  the  two  sintering  mechanisms  lead 
to  two  very  different  structural  evolution  signatures.  In  Figure  4,  we  show  the  average  pore 


%  TD 


Fig.  4.  /?efr(0)  as  a  function  of  percent  theoretical  density  for  silica 
and  alumina. 

radius  as  a  function  of  densification  for  the  two  systems.  The  alumina  results  are  familiar  from 
Figure  1 .  The  silica  results  show  that  the  average  pore  size  increases  and  then  decreases  in  the 
intermediate  stage  instead  of  remaining  constant.  This  is  because,  in  the  glassy  system,  the 
aggregates  first  densify,  removing  the  smaller  pores  inside  the  clusters.  This  leaves  behind 
larger  spaces  (pores)  between  the  clusters.  The  peak  radius  identifies  the  point  in  the 
densification  at  which  all  of  the  small  pores  have  been  removed.  After  that,  the  pore  size 
decreases  rapidly.  In  addition,  there  is  a  morphological  change  in  the  pore  shape  that  goes 
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from  interconnected  cylindrical  channels  to  interconnected  flat  (disklike)  pores  in  the  final 
stage.  The  transition  to  collapsed  porosity  in  silicia  li  currently  under  further  investigation. 
This  densification  behavior  is  attributed  to  the  lack  of  constraints  associated  with  the  absence 
of  a  grain  boundary  in  the  glassy  system.  The  density  of  the  green  body  plays  a  much  less 
important  role  when  viscous  flow  i*-  the  operative  sintering  mechanism. 

Further  research  on  processing/microstructure  relationships  will  include  silicon 
nitride  sintering  as  a  natural  progression  arising  from  the  program  just  completed.  The  new 
research  will  be  directed  toward  examining  conventional  silicon  nitride  processing  using  variety 
of  sintering  aids,  microwave-assisted  reaction-bonded  silicon  nitride  which  is  not  sintered  and 
involves  no  sintering  aids,  and  nanophase  processing  of  silicon  nitride  which  also  involves  no 
sintering  aids.  The  goal  will  be  to  understand  the  underlying  mechanisms  for  silicon  nitride 
fabrication  leading  to  optimized  reproducible  processing  of  this  promising  ceramic. 

4.  SYNTHESIS  AND  X-RAY  CHARACTERIZATION 

The  sol  gel  method  is  one  of  the  most  active  areas  of  ceramics  research.  Part 
of  the  synthesis  and  x-ray  program  supported  by  the  ARO  was  devoted  to  exploring  the  sol  gel 
method  with  regard  to  new  possibilities  in  the  preparation  of  alumina.  Earlier  investigations 
focussed  on  boehmite  (AfOOH),  its  conversion  to  a-alumina,  and  the  subsequent  sintering 
behavior.  In  this  program  we  have  synthesized  and  evaluated  a  hydrous  vitreous  alumina 
precursor  (HVAP).  These  precursors  have  been  prepared  both  with  and  without  the  addition 
of  a-alumina  seed  particles,  and  their  behavior  as  they  are  converted  to  corundum  (a-alumina) 
has  been  compared  to  that  of  microcrystalline  boehmite  precursors.  Our  x-ray  diffraction  and 
differential  thermal  analysis  data  indicated  that  the  HVAP  material  converts  to  corundum  at  a 
lower  temperature  than  the  boehmite.  Moreover,  me  presence  of  a-alumina  "seed"  particles 
in  these  precursors  lowers  the  respective  transition  temperatures  even  further.  A  comparison 
of  x-ray  diffraction  intensities  from  seeded  and  un„.:eded  samples  suggests  that  the  presence  of 
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seed  particles  has  an  effect  on  the  formations  of  -A£-OH-Af-  and  -Af-O-Af-  networks.  Low 
alumina  recoveries  from  gravimetric  alumina-content  assays  on  the  HVAP  indicate  that  it 
chemically  bonds  water  and  possibly  alcohol  into  its  structure.  In  contrast,  the  boehmite 
precursors  have  alumina  yields  more  consistent  with  theoretical  predictions.  We  conclude  that 
the  HVAP  materials  offer  a  new  route  to  synthesis  of  alumina,  and  further  research  is  planned 
to  continue  to  explore  these  possibilities. 

As  part  of  this  research  program,  we  have  measured  x-ray  diffraction  profiles 
during  the  in-situ  sintering  of  alumina.  (See  Figure  5.)  The  peaks  at  20  =  40°  and  58°  arise 
from  the  molybdenum  heating  strip.  The  sintering  process  causes  the  width  and  the  intensity 
of  the  sample  diffraction  profiles  to  change  as  the  temperature  increases.  The  initial  submicron 
size  of  the  crystallites  (Sumitomo,  HKP50,  0.2/im)  causes  the  diffraction  peaks  to  be  broad. 
During  the  early  stage  of  sintering,  crystallites  grow  out  of  this  range  and  the  diffraction 
profiles  become  narrower.  A  second  effect  is  extinction,  which  noticeably  reduces  the 
diffracted  intensity  of  the  strongest  lines  as  the  crystallites  coarsen  to  sizes  larger  than  those 
which  can  be  followed  by  peak  broadening.  Both  effects  are  sometimes  subtle  alterations  of 
the  diffraction  profiles.  The  analysis  of  this  data  requires  that  the  diffuse  thermal  scattering 
be  removed  and  the  effects  of  extinction  be  derived  from  the  fitting  of  the  profile  shapes.  Such 
analyses  are  underway  and  the  grain-size  evolution  data  is  expected  to  be  complementary  to  the 
pore-size  evolution  results  described  earlier. 

Finally,  we  describe  a  project  in  which  we  have  investigated  the  pore-size 
distribution  in  final-stage  sintering  by  means  of  small-angle  x-ray  scattering.  The  values  of 
the  median  size  and  the  width  of  the  pore-size  distributions  are  of  particular  interest  in  gaining 
an  understanding  of  the  processes  which  take  place  as  a  ceramic  body  approaches  full 
theoretical  density.  The  distribution  of  pore  sizes  is  related  to  distribution  of  dihedral  angles 
between  grains,  which  are  in  turn  directly  related  to  the  driving  force  for  sintering,  i.e.,  the 
surface  free  energy.  This  research  was  stimulated,  at  least  in  part,  by  our  neutron  scattering 
results  which  indicated  that  the  median  size  of  the  few  remaining  pores  in  final  stage  sintering 
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always  increased  as  density  increased.  We  measured  samples  of  alumina,  both  doped  with 
MgO  and  undoped,  that  were  between  97%  and  more  than  99%  dense.  The  high  resolution 
x-ray  results  showed  that  the  mean  pore  radius  continues  to  increase  as  density  increases, 
confirming  and  extending  the  neutron  results.  In  addition,  it  was  found  that  the  pore-size 
distribution  broadens  significantly  (Figure  6)  as  density  increases.  Scanning  electron 
micrographs  (Figure  7)  of  these  specimens  support  the  scattering  results  and  show  a 
concomitant  growth  of  the  alumina  grains.  Lastly,  for  a  given  sample  density,  the  average 
grain  and  pore  sizes  are  markedly  smaller,  and  their  size  distributions  are  narrower,  in  the 
MgO-doped  case.  This  study  gives  statistical  validity  to  earlier  research^  which  suggested  that 
the  distribution  of  dihedral  angles  between  grains,  and  thus  the  grain  and  pore  size 
distributions,  are  narrower  in  MgO-doped  alumina  than  undoped  alumina  at  —99%  dense. 


Fig.  6.  Volume  size  distributions  V(D),  of  the  pores  derived  using  the  maximum  entropy 
p^edure,  plotted  as  a  function  of  diameter,  D,  for  the  MgO-doped  alumina  samples  For 
clarity,  smoothed  lines  have  been  drawn  through  the  data  points. 
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Fif.  7.  Scanning  electron  micrographs  of  undoped  alumina  samples  13-7  (97.(Kf  TD).  13-8  TDl  and  of  MgO-dopcd  samples  D-6 

(97.5%  TD)  and  D-7  (99.0%  TD). 
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ABSTRACT 

Multiple  small  angle  neutron  scattering  (MSANS)  wae  used  to  follow 
tha  evolution  of  the  pore  size  distribution  In  a-Al^O^  through  the 

Intermediate  and  flnel  etegee  of  sintering.  This  technique  makee  It 
poeslble  to  determine  nlcrostructure  morphology  In  the  0.06  to  10  /ua  size 
regime  under  conditions  where  the  total  scattering  croes-sectlon  le 

dominated  by  elastic  Incoherent  multiple  ecatterlng,  as  Is  often  the  case 
for  ceramlce.  The  MSANS  results  Indicate  an  Initial  decrease  In  the 

effective  pore  radius  from  0.19  /im  at  57%  of  theoretical  denelty  (TC)  to 
0.17  ;im  at  79%  TO.  Ae  the  sample  density  further  Increased,  there  was  e 

transition  region  after  which  the  effective  pora  radius  grew  rapidly  to 

2:0.5  fim  et  98%  TD.  Standard  Porod  enalysle  on  ecatterlng  by  the  eame 
eamples  eleo  Indlcatee  e  transition  betwaan  the  Intermediate  end  final 
etegee  of  proceeelng.  Both  sets  of  results  support  a  topological  model  of 
e Inter Ing  In  which  the  Interconnected  pore  network  In  the  Intermediate  etege 
of  elnterlng  decaye  In  a  eteble  manner.  Tha  poree  become  fewer,  while 
retaining  the  same  diameter,  as  denslf Icatlon  proceeds.  In  the  final  etage, 
leoleted  pores  remain. 


INTROOUCTION 

Knowledge  of  the  micro structure  evolution  es  a  function  of  thermal 
proceeelng  la  Important  for  the  development  of  procees  modele  In  ceramlce. 
For  dilute  concentrations  of  scatterere  between  1  na  and  100  nm,  email  angle 
neutron  scattering  (SANS)  Is  Invaluable  for  the  determination  of  bulk  mlcro- 
etructural  parametere.  The  complementary  MSANS  technique,  first  deecrlbed 
by  Berk  end  Hardmen'Rhyne  le  applicable  when  the  scetterers  are  In 
tha  0.06  ;ia  to  10  >im  elze  range.  In  such  cases,  the  total  ecetterlng  crose- 
eactlon  le  dominated  by  elaetlc  Incoherent  multiple  ecatterlng.  Effective 
(volume 'Weighted  average)  radii  of  the  scatterere  era  extracted  by  meaeurlng 
the  wavelength  dependence  of  the  ecatterlng  curvae  et  low  ecatterlng  vector. 
Thus,  It  becomee  poeelble  to  Investigate  thicker  end  denser  materlale  than 
Is  poeslble  with  SANS.  Of  the  many  methode  evallable  for  the  determination 
of  mlcrostructura  morphology,  MSANS  and  SANS  have  the  adventege  of  being 
able  to  measure  both  open  end  closed  pores. 

In  this  work,  an  Inveetlgatlon  of  the  micro structure  evolution  of  m- 
Al^Oj  as  a  function  of  thermel  proceeslng  wae  conducted.  Thle  le  a  eystem 

for  which  the  procese  models  (3]  ere  comparetlvely  we 11 -under stood,  and  thue 
It  cen  eerve  ae  e  model  system  for  the  new  scetterlng  techniques.  An 
aarllar  etudy  of  mlcroetructura  evolution  ae  a  function  of  thermal 
proceeelng  In  porous  elllca  [4]  hae  demonstrated  that  MSANS  and  SANS 
effectively  cover  tha  full  range  of  relevant  mlcroetructura  elzee.  For  tha 
alumina  eyetem,  MSANS  wee  used  throughout  the  Intermediate  end  final  etegae 
because  tha  pora  elzae  were  lerga.  Standard  Porod  [5]  enalysls,  performed 
on  the  high  engle  portion  of  the  scattering  curves,  wae  used  to  determine 
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tocal  ecaccerlng  eurface  areae,  fron  which  effective  (surface  area-weighted 
average)  pore  radii  were  calculated. 


HSANS  THEORY 


The  Interecclon  of  neutrone  with  ■aCter  le  deCemlned  by  the  phaee 
ahlfc,  V,  chat  e  plane  wave  undergoee  In  traversing  a  particle  of  radlue,  R. 
Thle  phaee  ehlft  dcpende  on  the  Index  of  refraction  of  the  particle  or  void, 
An,  relative  to  the  ecatterlng  Mtrlx  euch  that  t'-(4e/l)hnR,  where  X  le  the 

2 

neutron  wavelength.  An-dbl  /2x,  where  Ab  le  the  relative  ecatterlng  length 
density,  or  contraet,  of  the  particle  or  void: 


Ab  - 


cell  I'  cell 


-  b 


Mtrlx' 


(1) 


Thus,  the  phaee  ehlft  w  le  given  by 

w  -  2AbRl,  (2) 

■eenlng  that  It  la  directly  dependent  on  the  Bacerlsl  contraet,  the 
dimensions  of  the  scatterer  and  the  neutron  wavelength. 

If  w  «  1,  the  measurement  le  In  the  conventional  SANS  diffraction 
regime.  The  form  of  the  ecatterlng  curve  as  e  function  of  the  ecatterlng 
wavevector  Q  (where  |Q|  -  Int/X,  and  t  le  the  ecatterlng  angle)  la 
Independent  of  neutron  wavelength  and  depends  only  on  the  particle 
dimensions.  On  the  other  hand.  If  i/  »  1,  the  experiment  Involves  multiple 
refraction,  where  every  particle  ecatters  and  geometrical  optics  applies. 

2 

The  scattering  curve  depends  on  neutron  wavelength  ae  1  .  When  i/  s  1,  the 
MSANS  diffraction  formalism  applies  Che  scattering  curve  broadens 

as  e  function  of  neutron  wavelength.  The  predicted  neutron  Intensity  le 

2 

approxlmetely  gausalen  a.  Q  -  0  and  hae  a  width  nearly  proportional  to  X  . 
Even  though  Che  wavelength  dependence  Is  the  eame  es  that  for  multiple 
refraction,  the  feet  that  w  5  1  allows  en  effective  particle  radius, 
R^^£(0),  to  ba  determined  from  the  wavelength  dependence  of  the  radius  of 

curvature,  r^,  at  Q  -  0  of  the  scattering  spectrum. 

A  Mesure  of  the  amount  of  multiple  scattering  Is  given  by  the 

parameter  z  which  le  the  sample  thickness,  z,  divided  by  the  statistical 
mean  free  path  length,  Jt. 

i  -  z/1  -  1.5  #z(AbA)^R^jj(0),  (3) 

where  ^  la  the  volume  fraction  of  acatterers  and  the  dimensions  of  the 

parameters  are  expressed  In  centimeters,  z  Is  related  to  the  scattering 
radius  of  curvature  via  an  empirical  relation  [2] 

r^tQRJ  -  0.926U  for  z  >  5.  (4) 


For  z  >  10,  the  egreeMnt  with  the  theory  Is  better  than  0.1%.  The  enalyela 
requlree  only  the  sample  thickness,  the  wavelength  dependence  of  the  radius 
of  curvature,  the  scattering  voIum  fraction  (which  Is  e  reflnable 
parameter)  end  the  scattering  length  density  to  derive  large 

Q,  the  scattering  Is  determined  by  the  single  scattering  Porod  lew  (5], 
which  applies  even  when  the  region  near  Q  -  0  la  dominated  by  multiple 
scattering  [1). 


18 


63 


ALUMIKA  SAKFUS 

Tha  suplae  vara  Saikowski  alunlna  powdar  (CR6},  sllpcasc  Co  an  avarega 
graan  daneity  of  33%  TD.  Tha  allp  conslsced  of  43.3  vol.%  eluntna  powder 
■ixad  into  wacer,  along  with  0.3  wc.%  of  a  polyalaccrolyca  dtspersanc 
(Darvan  C) .  Tha  suspanslon  was  ulCrasonicaCad  for  13  ain  and  allowad  Co  aac 
ovamighc  bafora  baing  poured  Inco  aoulda  wich  caflon  rings  on  Cop  of 
plascer-of-paris  blocks.  Afcar  33  Binuces,  Cha  aoulds  wars  ranovad  froa  che 
blocks  and  saalad  In  a  container  overnight.  Tha  parclally-drlad  saaplas 
ware  than  ranovad  froa  che  Mulds  and  aaaled  In  concalnars  for  2  days.  Tha 
concalnera  were  Chan  opened  co  tha  air  for  an  addlclonal  3>7  dsye  before  a 
aerlea  of  ovan*drylng  aCapa  ware  perfomad  ee  follows:  2  days  at  64*C,  1 
day  ac  81*C,  1  day  at  102‘C  and  1  day  ac  116*C.  afcar  which  a  dry,  eolld 
green  body  la  obtained.  In  chls  fashion,  a  faally  of  eight  nearly* Identical 
aattplas  were  prepared. 

During  sintering,  aaaplas  were  raaoved  one  by  one  froa  the  furnace 
co  yield  exaaples  of  aacarlal  between  36.3  end  98%  TD.  Figure  1  Is  a  time* 
Ceaperatura  curve  for  Cha  alntarlng.  Each  circle  rapraaante  a  eaapla 
raaoved  froa  che  furnace  and  Che  ntosber  alongside  Is  Che  aeasurad  %TD. 


NEUTRON  SCATTERING  MEASUREMENTS 

Maesureaants  ware  aada  at  che  20  MW  research  reactor  at  che  National 
Inscicuca  of  Standards  and  Technology  (NIST,  foraarly  Che  National  Bureau  of 
Standards).  The  NIST  SANS  facility  (6}  includes  a  velocity  aelactor  for 
chooaing  the  aean  Incident  wavelength.  X,  of  Cha  neutrons  with  a  spread, 
tX/X  •  0.23.  Copious  neutron  fluxes  can  be  obtained  in  Che  wavelength 
region  baewaan  0.3  and  2.0  na  due  Co  a  cryogenic  aodarator  installed  in  Che 
reactor  core.  Tha  neutron  beaa  is  colllaatad  wich  a  12  aa  aperture  which 
follows  tha  velocity  aalactor  and  a  8  aa  aperture  which  pracadea  tha  eaaple 


a. 

s 


uo  i«e  i«o  UD 
Tia*  (a<n) 


Figure  1.  Sintering  ciaa- 
taaparatura  curve  for  the 
eluaina  eaaplas  sintered  to 
densities  ranging  froa  56.3 
to  98%  TD. 
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poeltlon.  ScaCCered  ncucrone  ara  detected  by  e  6^  x  64  cm  posltlon- 
eeneitlve  detector  divided  Into  128  x  128  plxele.  The  data  are  circularly 
averaged  to  produce  one-dlmenelonal  intensity,  I(Q),  versus  Q  curvee.  For 
Porod  analyels,  ecatterlng  cross-eactlons  ara  placed  on  an  absolute  scale  by 
normalizing  to  that  of  e  standard  silica  gel  sample  measured  under  the  eame 
conditions . 


RESULTS  Ain)  DISCUSSION 

MSANS  maasurementa  ware  performed  on  each  sample  at  a  minimum  of  S 
different  neutron  wavalengthe.  The  scattering  curves  exhibited  the  expected 
broadening  with  Increasing  neutron  wavelength.  Causalan  flte  were  made  to 
each  curve  In  the  range  extending  from  90t  of  1(0),  the  intensity  at  Q-0,  to 
AS%  of  1(0) .  The  scattering  curvee  no  longer  fit  wall  to  a  gausalan  in  the 
higher  Q  region,  where  the  ecattaring  Intensity  falls  below  40%  of  1(0). 
For  each  fit,  a  radius  of  curvature  wee  extracted  from  the  etandard 

deviation,  r^  >  ■/i/2c.  Blank  epcctra  (with  no  acatterar  in  place) ,  maaeurad 

at  aach  wavelength,  ware  aleo  fitted  in  order  to  make  a  correction  for  the 
instrument  function.  The  meaeured  values  of  r^  (aymbola)  along  with  the 

theoretical  valuae  determined  using  the  HSANS  formellem  (solid  lines)  are 
plotted  vareue  neutron  wavelength  in  Figure  2  for  each  eampla.  The 
dansltlaa  meaeured  volumetrlcally  and  those  determined  by  the  neutron 
ecattaring  are  shown  in  the  inset  to  the  figure.  The  derived  values  of 
R  ££(0)  fit  the  theory  well  for  all  eamplca  meaeured. 


Figure  2.  Heeeursd  radii 
of  curvature  (aymbols)  end 
the  corresponding  theo¬ 
retical  radii  (eolld  llnae) 
as  a  function  of  neutron 
wavelength  for  all  samples 
meaeured.  The  eampla 
densities  meaeured  volum- 
strlcally  ara  shown  in  the 
inset,  along  with  those 
determined  by  tha  neutron 
scattering  meaeiirementa ,  If 
different. 


o.i  o.t  i.a  i.z  i.«  M 
Wavelength  (nm) 
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R^^^(O)  was  darlvad  from  Flyura  2  for  each  alumina  temple.  A  plot  of 

R^^^(O)  as  a  funcCion  of  ITD  le  illuatraced  in  Figure  3,  along  vlch  the 

evarege  pore  radius.  darlved  from  single  parclcla  Porod  measure- 

mente  assuming  spherical  scaccerers.  The  MSANS  decermlna clone,  which  are 
strongly  weighted  to  Che  larger  volume  particles  [A],  Indicate  chat  Che 
average  pora  radius  for  the  distribution  within  the  alumina  decreaaee 
slightly  from  0.19  pm  to  0.17  pm  In  the  Intermedlace  sintering  etaga  and 
Increases  sharply  to  >0.5  pm  In  the  final  acage.  Tha  Porod  de terminations , 
which  are  weighted  to  the  smaller  volume  parclclea,  Indicate  an  average 
pore  radius  chat  le  slowly  Increasing  during  Che  Intermediate  stage  of 
sintering  and  more  rapidly  Increasing  during  Che  final  stage.  Both  results 
Indicate  a  transition  in  the  mlcrostructure  morphology  between  85  and 
95*  TD.  The  quantitative  differences  arise  In  part  becausa 

darlvad  from  a  meaaure  of  the  volume  fraction  of  voids  In  the  system  and 

R^^£(m)  from  a  meaaure  of  the  total  eurfaca  area  In  the  spaces  (or  necks) 

between  tetrahedrally  positioned  particles.  Mercury  poroelmetry 
meaauramenta  (R.A.  Page,  unpubllehed  resulte)  performed  on  the  samples  below 
85t  TD.  where  open  porosity  Is  axpectad,  agree  qualitatively  with  both  sacs 
of  SANS  results  and  the  average  pore  radii  are  quantitatively  closer  to  the 
R^££(«>  values. 

The  Important  result  of  both  the  SANS  and  poroslmetry  messuremants  la 
Chat  the  effective  pora  radius  remslne  relatively  constant  through  the 
Intermediate  stage  of  sintering.  This  le  consistent  with  Intemedlate  etaga 
sintering  models  In  which  topological  decay  of  an  Interconnected  pora 
network  le  taking  place  [7].  Throughout  Che  Intermediate  stage  of 
sintering,  the  ratio  4/^^,  or  volume  fraction  to  surface  area  par  unit 

volume,  la  constant  [8],  Consequently,  Che  pore  channels  retain  a  constant 
dlametar  as  danslflcation  proceeds. 

Tha  affective  pore  radii  measured  with  MSANS  correspond  Co  the  cross* 
section  of  Che  pore  channels.  Tha  average  length  of  the  pora  channels  is 
>10  pm,  which  is  beyond  the  range  of  the  MSANS  technique.  Thus,  R  ^^(O) 


Perctnt  of  Theoretical  Density 

Figure  3.  Effactlva  pora  radii  for  tha  alumina  samples  as  a  function  of 
percent  thaoratlcal  density  dacarmlnad  from  the  MSANS  formalism 

and  from  the  single  parclcla  Forod  meaaurements  ». 
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rculna  ralatlvaly  constant  oa  danslflcatlon  proceeds  and  the  pore  channel 
natvorlc  dacaya. 

Before  the  channala  have  disappeared  altogether,  a  transition  into  the 
final  sintering  atage  occura  in  which  the  poroalty  beconas  isolated.  Pore 
iaolatlon  poaalbly  occura  due  to  a  rupture  or  ’pinching  off*  of  the  channala 
as  they  can  no  longer  mslntaln  their  length  while  at  the  saisa  tiae  raaalnlng 
atabla  against  the  increasing  aintarlng  force  [7].  The  Isolated  porea  are 
■ore  spherical,  having  dlaaaters  tdiich  are  greater  than  the  croas- sectional 
dlaisatara  of  the  original  channels,  yet  saall  enough  to  ba  within  the  MSAMS 
range.  Aa  the  nuabar  of  iaolatad  pores  Increaaaa,  tha  larger  diaanaion 
■ekes  an  increasing  contribution  to  the  total  acstterlng,  until  It  finally 
dominates  the  MSAMS  crosa>aectlon.  Thus,  Increases  gradually  as  tha 

transition  from  intaraadiata  to  final  sintering  atage  begins.  Late  In  tha 
final  stage,  a  more  rapid  increase  is  seen  aa  the  aicrostructure  becoaea 
doainated  by  the  Isolated  porea. 


CONCLUSIONS 

The  transition  from  one  sintering  stage  to  another  in  sintered  alualna 
was  detectad  by  saall  sngle  neutron  scattering.  The  affective  radius  of  tha 
pore  size  distribution,  aa  neasured  by  both  ■ultlple  ssall  angle  neutron 
acstterlng  and  single  particle  Porod  scattering,  raaalned  fairly  constant 
throughout  Che  InCcracdlaca  scsge  of  aintarlng.  with  a  rapid  Incraaaa  aa 
densif Icatlon  proceeds  beyond  85%  of  thaorecleal  density.  The  results  are 
consistent  with  sintering  nodals  which  aquata  the  InCeraadiate  stage  with 
tha  decay  of  the  topological  network  of  chsnnals  and  Junction  poras  in  tha 
aystaa.  Before  coaplaCa  decay  can  occur,  a  Cranalcion  into  the  final 
sintering  atsga  occurs  in  which  the  porosity  becoaea  iaolatad  spherical 
voids. 
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Smatl-angle  X-ray  scattering  was  used  to  follow  the  evolu¬ 
tion  of  the  pore  size  distribution  during  final-stage  sintering 
of  alumina  and  of  alumina  doped  with  0.25  wt%  magnesia. 
Tbe  volume-weigbted  (Guinier)  results  indicate  that  the  ef¬ 
fective  size  of  the  largest  pores  increases  as  the  body  goes 
from  97%  to  more  than  99%  dense.  The  surface-area- 
weighted  (Porod)  results  show  that  the  median  size  of  the 
smallest  pores  decreases  slightly  over  the  same  density 
range.  Taken  together,  these  data  indicate  that  the  pore  size 
distribution  becomes  broader  as  final-stage  densification 
proceeds.  This  was  confirmed  by  a  maximum  entropy  analy¬ 
sis,  which  was  used  to  derive  pore  size  distributions  directly 
from  the  data.  Finally,  the  evolution  of  the  pore  size  distri¬ 
butions  io  alumina,  with  and  without  sintering  aid,  were 
compared.  (Key  words:  alumina,  sintering,  pore  size  distri¬ 
bution,  X-ray.  magnesia.) 

I.  Introduction 

The  microstructurc  evolution  of  the  porous  phase  during 
final-stage  sintering  of  alumina  and  of  alumina  with 
0.25  wt%  MgO  sintering  aid  has  been  studied  by  means  of 
high- resolution  small-angle  X-ray  scattering  (SAXS).  The 
values  of  the  median  and  the  width  of  the  pore  size  distribu¬ 
tions  arc  of  special  interest  in  gaining  an  understanding  of  the 
processess  which  take  place  as  the  ceramic  body  approaches 
full  theoretical  density.  In  particular,  the  distribution  of  pore 
sizes  is  related  to  the  distribution  of  the  dihedral  angles  be¬ 
tween  grains,  which  are  in  turn  directly  related  to  the  driving 
force  for  sintering,  i.e..  the  surface  free  energy. 

Earlier  studies' of  the  microstructure  evolution  of  alumina 
as  a  function  of  sintering  were  performed  using  multiple 
small-angle  neutron  scattering.  In  that  work,  both  the  inter¬ 
mediate  and  final  stages  of  sintering  were  studied  and  the 
results  indicated  that  the  effective  pore  radius  is  constant 
during  the  intermediate  stage,  independent  of  the  presence  of 
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MgO.  During  that  stage,  the  value  of  the  effective  pore  radius 
was  shown  to  depend  on  the  green  body  density  or,  equiva¬ 
lently,  on  the  connectivity  in  the  compact  before  sintering. 
When  final-stage  sintering  begins,  the  open  interconnected 
porosity  of  the  intermediate  stage  breaks  up  into  isolated 
porosity,  and  the  effeaive  pore  size  increases  abruptly.*  It 
was  found  that  the  role  of  MgO  as  a  sintering  aid  lies,  at  least 
in  part,  in  prolonging  the  stability  of  intermediate-stage  sin¬ 
tering  such  that  the  body  achieves  greater  density  before  the 
transition  to  isolated  porosity  takes  place,  leading  to  a  finer- 
grained  finished  ceramic  in  the  case  of  the  densified  material 
containing  MgO.  Contrary  to  expectation,  the  measured  ef¬ 
fective  pore  radius  was  never  observed  to  decrease  as  densiiy 
increased,  even  for  material  that  had  reached  98%  of  theoreti¬ 
cal  density  (TD). 

The  present  study  was  undertaken  to  measure  the  evolution 
of  the  pore  size  distributions  during  late  final-stage  densifica¬ 
tion.  This  kind  of  information,  which  can  be  derived  from 
single-particle  scattering,  is  not  generally  available  from  multi¬ 
ple  scattering  measurements  unless  the  shape  of  the  distribu¬ 
tion  is  known.  A  high-resolution  SAXS  camera.’*  which  can 
measure  sizes  between  3  nm  and  I  ^m,  was  used  to  measure 
samples  of  alumina  and  alumina  doped  with  0.25  wt%  MgO 
with  densities  ranging  from  97%  to  more  than  99%  of  theo¬ 
retical.  The  effective  mean  pore  radius  was  found  to  increase 
as  density  increased,  confirming  the  earlier  multiple  neutron 
scattering  results.  In  addition,  it  was  found  that  the  pore  size 
distribution  broadens  significantly  as  density  increases.  This 
behavior  can  be  compared  to  the  trends  that  are  observed’  in 
the  distributions  of  dihedral  angles  between  grains. 

11.  Samples 

The  samples  (13-series  and  8-series)  were  made  from 
Baikowski  alumina  powder  (CR6,  6  m’-g'')  and  slip  cast  to 
green  densities  of  approximately  53%  of  theoretical,  where 
the  theoretical  density  of  alumina  is  3.985  g  cm'L  The  slip 
consisted  of  43.5  vol%  alumina  powder  mixed  into  water, 
along  with  0.5  wt%  of  a  polyelectrolyte  dispersant.'  The  sus¬ 
pension  was  agitated  ultrasonically  for  15  min  and  allowed  to 
settle  overnight  before  being  poured  into  molds  with  Teflon’ 
rings  on  top  of  plastcr-of-paris  blocks.  The  following  efforts 
were  made  to  guarantee  that  the  material  was  dried  slowly  to 
prevent  cracking.  After  55  min,  the  molds  were  removed  from 
the  blocks  and  sealed  in  a  container  overnight.  The  partially 
dried  samples  were  then  removed  from  the  molds  and  sealed 
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into  containers  for  2  days,  after  which  the  containers  were 
opened  to  the  air  fc  5  to  7  days.  Finally,  ••  series  of  oven- 
drying  steps  were  performed;  2  days  at  64°C,  I  day  at  SrC. 
I  day  at  102°C,  and  I  day  at  116°C.  At  this  point  a  dry,  solid 
green  body  was  achieved.  Compacts  (D-series)  were  also  pre 
pared  by  cold-pressing  from  the  same  powder  doped  with 
0.25  wt%  MgO.  Green  densities  of  approximately  of  the¬ 
oretical  were  obtained  for  these  doped  samples.  Doping  was 
achieved  by  adding  a  solution  of  MgfNOsh  in  distilled  water 
and  drying  overnight  at  200°C.  Earlier  neutron  scattering 
studies'  showed  that  the  green  body  density,  rather  than  the 
preparation  method,  determines  the  sample  characteristics. 

The  samples  were  heated  to  1600°C,‘”  at  which  time  90% 
TD  was  achieved.  With  the  furnace  temperature  held  con¬ 
stant  at  IbOO'C,  samples  were  removed  one  at  a  time  to  yield 
examples  of  sintered  microstructure  between  97%  and  more 
than  99%  TD.  The  characteristics  of  the  three  suites  of 
samples  are  given  in  Table  I.  Pore  volume  fractions  were  de¬ 
termined  volumetrically  with  an  uncertainty  of  0.3%.  Samples 
13-7,  13-g,  and  D-6  were  taken  from  samples  measured  earlier 
by  neutron  scattering.'-*  The  D-7A  sample  was  obtained  by 
further  sintering  a  piece  of  the  D-7  sample  to  99.7%  TD.  The 
samples  for  the  X-ray  studies  were  cut  out  of  the  larger  sam¬ 
ples  using  a  standard  3-mm-diameter  electron  microscopy  die 
and  then  polished  to  thicknesses  between  150  and  200  /xm. 

HI.  Experimental  Procedure 
(1)  SAXS  Instrument  and  Measurements 

The  small-angle  X-ray  scattering  measurements  were  made 
by  means  of  a  modified  Bonse-Hart-type  double-crystal  dif¬ 
fractometer'  on  the  X23Aj  NIST  Materials  Science  beamline 
at  the  National  Synchrotron  Light  Source,  Brook  haven  Na¬ 
tional  Laboratory.  The  double-crystal  instrument  offers  high 
angular  resolution  at  the  smallest  scattering  angles  where  the 
scattering  from  large  microstructure  is  to  be  found,  and  it 
covers  the  entire  relevant  size  range  for  the  alumina  samples. 

A  diagram  of  the  experiment  is  shown  in  Fig.  1.  The  energy 
of  the  incident  photons  is  chosen  using  a  monochromator  con¬ 
sisting  of  two  Geflll)  crystals.  The  size  of  the  X-ray  beam  is 
defined  by  slits  at  the  monochromator  exit.  Since  the  X-ray 
beam  is  nearly  parallel,  the  slits  are  adjusted  to  maximize  the 
coverage  on  each  sample.  The  incident  photon  intensity  is 
monitored  by  means  of  a  combination  ionization-chamber/ 
helium-transport  tube  after  the  slits  and  before  the  sample 
position.  Samples  were  placed  on  kapton  tape  and  suspended 
from  rings  which  were  supported  at  the  sample  position  by 
meaits  of  two  fixed  pins.  The  scattered  photons  were  ana¬ 
lyzed  by  two  additional  Ge(lll)  crystals,  mounted  so  that  they 
rotate  as  a  monolith  about  the  center  of  the  surface  of  the 
first  crystal,  and  were  detected  by  a  5.5-mm-diameter  photo¬ 
diode  detector.* 

Small-angle  scattering  measurements  were  performed  on 
the  alumina  samples  using  11-keV  photons.  Each  sample 
was  aligned  individually  in  the  X-ray  beam  and  the  beam 


definition  slits  were  adjusted  to  maximize  the  coverage  in 
each  case.  After  the  small-angle  scattering  from  the  sample 
was  measured,  the  empty  beam  (with  no  sample  in  place) 
was  measured  before  aligning  the  next  sample.  Thus,  each 
sample  measurement  is  paired  with  a  correspotiding  empty 
beam  measurement. 

(2)  SAXS  Data  and  Analysis 

The  raw  data  of  photon  counts  vs  analyzer  crystal  angle 
were  reduced  to  SAXS  intensity  vs  scattering  vector  q  - 
4ir  sin  where  6  is  one-half  of  the  scattering  angle  and 
A  is  the  photon  wavelength.  The  intensity  was  placed  on  an 
absc'ute  scale  by  a  primary  calibration  technique.*  Since  the 
double-crystal  instrument  involves  slit  collimation,  it  was  nec¬ 
essary  to  perform  a  geometrical  correction  which  was  made 
using  the  desmearing  algorithm  of  Lakc.^  The  desmeared  in¬ 
tensities,  I(q),  for  each  of  the  three  suites  of  samples  are  plot¬ 
ted  on  a  In  ll{q))  v%q‘  scale  in  Fig.  2.  Note  that,  in  each  'use, 
the  scattering  from  the  most  dense  samples  has  the  steepest 
slope  at  small  q  values  while  the  least  dense  samples  show  the 
least  curvature  in  the  same  q  range.  Such  curvature  is  charac¬ 
teristic  of  samples  which  contain  a  distribution  of  sizes  of 
scatterers.  The  slope  of  the  curves  at  small  q  is  determined  by 
the  size  of  the  largest  pores  in  the  size  durribution.  Since  the 
steepest  slopes  are  associated  with  the  most  dense  samples, 
the  maximum  pore  size  continues  to  increase  as  the  material 
approaches  full  theoretical  density. 

Since  the  alumina  samples  do  not  contain  a  monodisperse 
pore  size,  a  single  radius  of  gyration  could  not  be  derived 
from  a  conventional  Guinier*  analysis  in  the  small-<j  region. 
However,  for  the  purpose  of  comparing  samples,  fits  to  the 
desmeared  data  were  made  in  the  same  limited  range 
5  X  10*'  £  s  6.5  X  10"'  nm  ‘  to  the  Guinier  equation 

/(<?)  =  /(O)exp(-(i- /?, 73)  (1) 

where  /I,  is  the  radius  of  gyration  of  the  larger  particles  and 
AO)  is  the  extrapolated  scattering  intensity  at  -  0  from  the 
range  of  scattering  vector  values  above. 

The  desmeared  scattered  intensity  beyond  q  -  0.03  nm  ' 
follows  the  standard  Porod  law'*  which  can  be  writter 

nq)  ^  Pq-'  +  B  (2) 

where  P  is  Porotfs  constant  and  B  is  a  background  term.  Both 
P  and  B  arc  obtained  from  linear  least -squares  fits  to  Iq*  vs  q* 
plots  for  each  sample.  The  total  scattering  surface  area  per 
unit  volume.  S/K  is  determined  directly  from  P  since 

P  =  2Tr{SSp)-S;V  (3) 

The  contrast,  Ap.  is  the  difference  in  X-ray  scattering  length 
density  between  the  particle  matrix  and  the  pores.  F^r  the 
alumina  samples,  Ap  is  equal  to  32.6  x  10'"  cm''.  An  average 
pore  radius,  Bp,  can  be  calculated  from  StV  if  the  pore  shape 
is  known.  After  the  background  term  is  subtracted  from  l(q), 
the  scattered  intensity  curves,  plotted  on  a  log  (A^))  ^s  log  (q) 
scale,  exhibit  a  slope  of  -4.  A  Porod  plot  from  sample  13-g, 


Table  I.  Alymina  Samples 


Simple 

Pore  volume 
fraction 

Green 

dcnsilv 

(■STD) 

Time 

ai  IMMTC  (mint 

Sample 

density 

(%TD) 

Sampie 

density 

(g-cni-  ’) 

Undoped  alumina 

13-7 

3.0 

53.0 

30 

97.0 

3.865 

138 

1.0 

52.9 

300 

99.0 

3.945 

8-10 

2.3 

52.8 

30 

97.7 

3.893 

8-8 

1.4 

52.9 

120 

98.6 

3  930 

8-14 

0.7 

52.9 

300 

99.3 

3.9.57 

M^>doped  Alumina 

2.5 

50.7 

7 

97.5 

3.885 

D-7 

1.0 

510 

3C 

99.0 

3.945 

D-7A 

0.3 

51.0 

180 

99.7 

3.973 

*Tht  uncertainty  on  il(  pore  volume  frKHon  valuei  is  ±0.3 
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Fig.  1.  High -resolution  smali-angie  X-ray  scattering  instrument. 


which  is  representative  of  Porod  scattering  from  all  of  the 
samples,  is  shown  in  Fig.  3.  The  corresponding  fitteo  slope  ts 
-4.02  i  0.06  in  the  range  0  03  s  <  0.11  nm  ‘. 

The  radius  of  gyration  obtained  from  Eq.  (1)  is  a  volume- 
derived  parameter  and  thus  is  weighted  to  the  larger  pores  in 
the  size  distribution  while  the  Porod  radius  is  a  surface-area- 
derived  parameter  and  is  weighted  toward  the  smaller  pores. 


Thus,  the  two  sets  of  results  taken  together  arc  a  qualitative 
measure  of  the  width  of  the  pore  size  distribution  for  each 
sample.  To  obtain  a  quantitative  description  of  the  distribu¬ 
tion  of  pore  sizes,  the  desmeared  SAXS  intensities  were  ana¬ 
lyzed  for  a  volume-weighted  distribution  (Ap)V(D),  where 
V[D)  dD  is  the  volume  fraction  of  pores  with  diameters  be¬ 
tween  D  and  D  +  dD.  The  size  distributions  were  derived 


(A)  (B) 


(C) 

rjg.2.  In  (/)  siq'  (Guinier)  plots  of  desmeared  SAXS  intensities  for  the  undoped  alumina  samples:  (A)  13- series  and  (B)  8-scries.  and  the 
MgO-doped  alumina  samples.  (C)  D-series. 
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log(q  nm'*) 


Fig.  3.  log  (/)  vs  log  (q)  (Porod)  plot  of  ihe  desmeared  SAXS  in¬ 
tensity  for  sample  13-8,  '^ich  is  a  typical  exam^e  of  Porod  scatter¬ 
ing  from  each  of  the  alumina  samples  studied.  The  solid  line  has  a 
fitted  slope  of  -4.02  z  0.06. 

using  a  maximum  entropy  technique'"  for  spherical  scatterers 
implemented  in  the  program  MAXE  obtained  from  the 
UKAEA  Harwell  Laboratory. 

The  maximum  entropy  algorithm  has  been  successfully  ap¬ 
plied  to  a  real  system  in  the  study  of  particle  size  distributions 
as  a  function  of  starting  chemistry  in  precursor  porous  silica 
samples."  In  this  earlier  work,"  it  had  been  emphasized  that 
it  is  important  to  apply  the  maximum  entropy  method  to  scat¬ 
tering  data  encompassing  a  q  range  sufficient  to  cover  the 
entire  size  range  of  the  pores  in  the  sample.  Otherwise,  the 
resultant  curve  fit  will  not  accurately  represent  the  true  pore 
size  distribution.  The  prior  assumption  used  in  the  maximum 
entropy  method  is  that  all  particle  dimensions  are  equally 
probable.  HeiKe  all  features  in  the  calculated  distributions 
are  only  produced  by  corresponding  statistically  significant 
features  in  the  scattering  data. 

One  measure  of  the  integrity  of  the  data  inversion  is  a 
comparison  of  the  pore  volume  fraction,  4>,  calculated  from 
maximum  entropy  with  the  volume  fraction  determined  volu- 
metrically.  Another  measure  of  the  integrity  of  the  data  is  a 
calculation  of  the  invariant'*  of  the  scattered  intensity.  The 
invariant,  often  represented  as  Q,  can  be  written  as 

Q  =  f  Hq)q^  dq  (4) 

When  l{q)  is  measured  on  an  absolute  scale,  Eq.  (4)  can  be 
rewritten  in  terms  of  the  scattering  volume  fraction  such  that 

Q  =  (5) 

Note  that  the  shape  of  the  pore  size  distribution,  or  equiva¬ 
lently,  the  shape  of  the  scattered  intensity  curve,  does  not 
enter  into  the  calculation  of  Q.  Whatever  the  shape  of  l{q)  vs 
q  for  a  given  value  of  the  value  of  the  invariant  does  not 
change,  provided  that  the  data  arc  sufficient  to  encompass  all 
pore  sizes  in  the  distribution.  Thus,  by  integrating  l{q}  q^ 
over  q,  the  volume  fraction  can  be  calculated  directly  from 
the  data  and  compared  to  that  measured  volumetrically. 

(9  SEM  Sample  Preparation  and  Measurements 

In  preparation  for  scanning  electron  microscopy  (SEM) 
measurements,  the  surfaces  of  the  specimens  were  ground 
with  successively  finer  SiC  grit  and  given  a  final  polbh  with 
1-nm  diamond  paste.  To  improve  the  visibility  of  surface  fea¬ 
tures,  the  samfAes  were  thermally  etched  by  heating  in  air  at 
ISOirC  for  IS  min  and  then  removed  at  temperature  and  al¬ 
lowed  to  cool  in  air.  The  samples  were  viewed  in  the  scan¬ 
ning  electron  microscope  at  both  12-  and  24-kV  accelerating 


Table  II.  SAXS  Guinier  Results* 


Sample 

Pore  volume 
fraction  (^) 

ft,  (M"’)' 

Undoped  alumina 

13-7 

3.0 

0.30  z;  0.02 

13-8 

1.0 

0.37  ±  0.02 

8-10 

2.3 

0.26  z:  0.01 

8-8 

1.4 

0.30  r  0.02 

8-14 

0.7 

0.31  z:  0.03 

Doped  alumina 

D-6 

2.5 

0.31  ±  0.01 

D-7 

1.0 

0.30  *  0.01 

D-7A 

0.3 

0.39  ±  0.02 

'Guinier  fits  were  always  made  in  (he  range  5  x  10' ' s  ^  s  6.J  x 
10*'  nm'*.  ’Radius  of  gyration. 


potential  and  at  0®  and  25®  till.  Magnification  was  checked 
with  NIST  SRM  4«4  (s/n  JY-43-cc). 

IV.  Experimental  Results 
(1)  SAXS  Results 

Values  for  the  radius  of  gyration,  R^,  obtained  from  a 
Guinier  analysis  in  the  data  range  5  x  10'^  S  ^  £  6.5  x 
10"’  nm'*  are  shown  in  Table  II.  These  results,  while  weighted 
toward  the  larger  pores  in  the  size  distribution,  are  not  taken 
to  represent  the  maximum  pore  size  since  the  curvature  of 
the  data  at  the  smallest  q  values  measured  is  still  increasing 
(see  Fig.  2).  which  indicates  that  still  larger  pore  sizes  may  be 
present.  The  Guinier  analysis  is  generally  valid®  when  the  q 
range  used  :s  centered  about  qRf  £  1.2.  In  this  case. 
qRf  —  1.5  at  the  smallest  q  value  measured.  Therefore,  the  q 
range  accessible  with  the  high-resolution  SAXS  instrument, 
while  adequate  to  detect  the  largest  pore  sizes,  is  not  suffi¬ 
cient  to  account  quantitatively  for  all  of  the  large  pores  in  the 
samples.  However,  these  results  can  be  used  to  place  a  lower 
limit  on  the  maximum  pore  sizes  present  in  the  samples  and, 
since  they  were  obtained  in  the  same  q  range  for  each  sample, 
they  can  be  meaningfully  compared. 

The  Rf  values  in  Table  11  show  an  increase  in  pore  size  with 
increasing  sample  density,  independent  of  the  presence  of 
MgO  sintering  aid.  This  is  illustrated  when  the  results  for  all 
three  suites  of  samples  are  plotted  together,  as  shown  in 
Fig.  4.  The  trend  observed  here  confirms  and  extends  the 


FIf.4.  Ptoi  of  pore  R,  (radius  of  gyration)  values  for  all  samples, 
obtained  from  a  Guinier  analysis,  as  a  function  of  sample  density. 
In  all  cases,  fils  were  made  in  the  data  range  i  x  10'^  S 
q  S  6.5  X  tO"^  nm"'. 
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Tabic  III. 

SAXS  Porod  Results 

Pore  volume 

qm,. 

P  t  *  |(l>  cm 

i/I  lx  III’ cm  l' 

Sample 

Iracuon  (^1* 

(fim  '1 

Inm 

Pr  lMmi‘ 

Undoped  alumina 

13-7 

3.0 

0.03 

0.11 

31.49  i  0.04 

4.716  X  0  006 

0.191  X  0.019 

13-8 

1.0 

0.03 

0.11 

11.67  r  0.02 

1.748  i  0.003 

0.172  X  0051 

8-10 

2.3 

0.03 

0.11 

25.38  =  0.02 

3.800  X  0.003 

0.182  X  0.023 

8-8 

1.4 

0.03 

0.11 

12.05  c:  0.02 

1.719  X  0.003 

0.233  X  0.053 

8-14 

0.7 

0.03 

0.11 

8.38  -  0.01 

1.254  X  0.001 

0.167  X  0.072 

Doped  alumina 

D-6 

2.5 

0.03 

0.13 

25.87  r  0.03 

3.874  X  0.005 

0.194  X  0.023 

D-7 

1.0 

0.03 

0.13 

11.06  X  0.02 

1.656  X  0.003 

0.180  X  0.055 

D-7A 

0.3 

0.03 

0.13 

4.89  X  0.02 

0.732  X  0.003 

0.123  X  0.123 

‘The  unceriiiiniy  on  nil  pore  volume  fraction  values  is  Porod  constant.  ‘Normaliaed  surface  area.  ^Porod  radius.  Spherical  scalterers  are  assumed. 


trends  observed  in  earlier  neutron  scattering  studies'"  where 
the  effective  pore  radius  was  observed  to  increase  monotoni* 
cally  with  increasing  sample  density  during  final-stage  sinter¬ 
ing.  In  addition,  it  was  shown  in  one  of  these  studies'  that  the 
effective  radius  of  pores  in  the  doped  sample  was  consistently 
smaller  than  that  of  pores  in  undoped  samples  at  the  same 
density.  Quantitatively,  the  pore  sizes  in  the  neutron  scatter¬ 
ing  studies  are  somewhat  larger  than  reported  here.  This  oc¬ 
curs  because  the  pore  volume  fractions  used  in  the  neutron 
scattering  studies  were  derived  directly  from  the  data  fitting 
procedure  and  are  somewhat  lower  than  the  volumetncally 
determined  values  used  in  this  work.  In  addition,  the  multiple 
small-angle  neutron  seattering  technique  accounts  for  the 
largest  pores  in  the  sample,  which  serves  to  increase  the  effec¬ 
tive  pore  radius  derived  from  those  studies.'  * 

Values  for  the  pore  radius  obtained  from  a  Porod  analysis 
of  the  SAXS  data  and  assuming  a  spherical  pore  shape  are 
shown  in  Table  III,  along  with  the  derived  Porod  constant,  P, 
and  the  normalized  surface  area  per  unit  volume,  S/V.  The 
largest  possible  linear  region  in  the  log  (/)  vs  log  (9)  curve 
beyond  q  =  0.03  nm"'  was  fit  in  each  case,  as  given  in  the 
and  columns  of  Table  III.  Figure  S  shows  the  PProd 
results  for  all  three  suites  of  samples.  The  large  error  bars  for 
the  samples  at  the  highest  densities  are  a  direct  result  of  the 
0.3%  uncertainly  on  the  pore  volume  fraction  values.  Con¬ 
trary  to  the  behavior  of  the  f?,  values,  the  values  decrease 
slightly  with  increasing  sample  density,  again  independent  of 
the  presence  of  MgO  in  the  sample.  Since  the  A,  values  are 
weighted  to  the  larger  pore  sizes  in  the  distribution  and  the 


Plot  of  pore  (Porod  radius)  values  for  all  samples,  ob¬ 
tained  from  a  Porod  analysis  and  assuming  a  spherical  pore  shape, 
as  a  function  of  samples  oensity.  Fits  were  made  in  the  largest  pos¬ 
sible  linear  daia  range  beyond  q  •  0.03  nm*'. 


Rp  values  are  weighted  to  the  smaller  pore  sizes,  the  behavior 
of  the  radii  in  Figs.  4  and  5  indicates  that  the  pore  size  distri¬ 
bution  broadens  with  increasing  sample  density. 

The  volume  size  distributions  f'(D),  derived  for  each 
sample  using  the  ma.ximum  entropy  algorithm,  are  plotted 
versus  diameter,  D.  for  each  suite  of  samples  in  Fig.  6.  The 
sample  densities  which  have  been  used  to  identify  the  curves 
are  those  which  had  been  determined  volumcirically.  The 
large  peaks  in  the  size  distributions  represent  the  centroids  of 
the  distributions,  while  the  smaller  ripples  may  be  due  to 
mathematical  artifacts.”*”  Therefore,  smoothed  lines  have 
been  drawn  through  the  points  for  clarity.  The  volume- 
weighted  mean  pore  radius,  Ry,  calculated  from  the  resultant 
size  distribution,  is  shown  for  each  sample  in  Table  IV.  These 
values  agree  quantitatively  with  the  average  pore  radii  deter¬ 
mined  from  the  earlier  neutron  scattering  studies,' *  provided 
that  the  volumetrically  determined  values  for  the  pore  volume 
fraction  are  used  in  all  cases. 

The  broadening  of  the  size  distributions  as  a  function  of 
densification  can  be  seen  in  Fig.  6  and  the  increase  in  the 
mean  pore  radius  can  be  seen  in  Table  IV.  This  table  also 
includes  the  values  of  the  standard  deviation  of  the  mean 
pore  radius  which  is  seen  to  increase  as  sample  density  in¬ 
creases. 

(2)  SEM  Results 

Scanning  electron  micrographs  for  undoped  samples  13-7 
and  13-8,  and  doped  samples  D-6  and  D-7,  are  presented  in 
Fig.  7.  The  grain  boundaries  and  pores  ean  clearly  be  seen. 
While  some  of  the  pores  are  larger  than  1  ^m,  most  are 
smaller.  In  addition,  most  of  the  pores  reside  on  grain 
boundaries  while  others  are  trapped  inside  large  grains.  The 
average  grain  size  in  the  undoped  samples  (13-7  and  13-8)  in¬ 
creases  dramatically,  along  with  an  increase  in  pore  size,  as 
the  sample  goes  from  97%  to  >99%  TD.  A  similar  inerease  in 
grain  and  pore  size  is  seen  in  the  doped  samples  (D-6  and 
D-7),  A  comparison  of  the  micrographs  of  Fig.  7  with  the 
pore  size  distributions  in  Fig.  6  shov/s  that  the  SEM  results 
are  consistent  with  the  SAXS  results  derived  using  the  maxi¬ 
mum  entropy  technique.  Note  that  the  region  of  sensitivity 
for  the  SEM  studies  differs  from  that  of  the  SAXS  studies. 
The  SEM  technique  measures  surface  features.  Therefore,  it 
is  well-suited  for  studying  the  largest  pores  in  the  samples 
since  the  necessary  polishing  and  etching  processes  tend  to 
remove  the  small  pores  from  the  sample  surfaces.  On  the 
other  hand,  the  SAXS  technique  measures  bulk  samples  so  it 
is  sensitive  to  the  smaller  pores  as  well.  Figure  6  shows  that 
the  majority  of  t'^e  pores  in  the  samples  fall  below  the  size 
range  readily  accessible  with  SEM. 

The  most  striking  feature  in  the  eomparison  between  the 
doped  samples  and  the  undoped  samples  is  that  the  average 
grain  size  in  the  doped  samples  is  markedly  smaller  than  that 
of  grains  in  the  undoped  samples  of  approximately  the  same 
density.  Since  pore  growth  is  associated  with  grain  growth, 
this  would  imply  that  the  average  pore  size  is  also  smaller  in 
the  latest  stages  of  sintering  for  the  MgO-doped  samples. 
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Diameter  (nin) 

(C) 

Fig.  <.  Volume  size  disiribulions  ViD).  derived  using  I  he  maximum  eniropy  procedure,  plotted  as  a  function  of  diameter.  D,  for  the 
undoped  alumina  samples:  (A)  13-series  and  (B)  8-series,  and  the  MgO-doped  alumina  samples,  (C)  D-series.  For  clarity,  smoothed  lines 
have  Dcen  drawn  through  the  data  points. 


V.  Discussion 

The  motivation  for  studying  the  pore  evolution  during  the 
late  stages  of  sintering  of  undoped  and  doped  alumina  with 
high-resolution  SAXS  was  the  result  of  the  earlier  neutron 
scattering  study'  of  the  same  material  which  suggested  that 


Table  IV.  SAXS  Maximam  Entropy  Results 


Sample 

Pore  volume 
CrKtion  (%) 

R.  (um)* 

Undoped  alumina 
13-7 

3.0 

0.245  (0.113) 

13-8 

1.0 

0.336  (0.172) 

8-10 

23 

0.225  (0.105) 

SS 

1.4 

0.236  (0.114) 

8-14 

0.7 

0.280  (0.161) 

D(g^  alumina 

23 

0.225  (0.096) 

D-7 

1.0 

0.272  (0.119) 

D-7A 

0.3 

0.292  (0.130) 

*VW«iiK-wcigft(*d  mean  radius.  Sraodard  deviackias  of  the  mean  reditu 
are  given  ia  ptrenlbeaes. 


the  effective  pore  radius  continues  to  increase  with  increasii^ 
sample  density,  even  when  densities  reach  98%  TD.  The  neu¬ 
tron  scattering  results  for  the  undoped  13-  and  doped  D-series 
samples  are  reproduced  in  Fig.  E,  where  the  effective  pore 
radius  is  plotted  versus  %TD.  This  figure  shows  that  the  mean 
pore  radius  during  intermediate-stage  sintering  is  determined 
by  the  green  body  density  or,  in  other  words,  by  the  original 
connectivity  in  the  green  state.  The  transition  f^rom  interme¬ 
diate  stage  (where  the  porosity  consists  mainly  of  intercon¬ 
nected  channels)  to  final  stage  (where  the  porosity  becomes 
isolated)  apparently  occurs  between  86%  and  92%  TD,  where 
the  effeaive  radius  begins  to  increase  rapidly  with  increasing 
sample  density.  This  transition  occurs  at  a  greater  density  in 
the  case  of  the  MgO-doped  sample,  suggesting  that  the  rote  of 
MgO  as  a  sintering  aid  lies,  at  least  in  part,  in  extending  the 
stable  intermediate  sintering  stage,  thus  allowing  the  body  to 
achieve  greater  density  before  final-stage  sintering  b^ins. 

Altho^h  there  are  large  uncertainties  in  the  radii  of  the 
denser  final-stage  samples,  the  effective  radius  of  the  un¬ 
doped  alumina  sam{^e  is  consistently  larger  than  that  for  the 
doped  sample  at  the  same  density.  This  effect  might  be  even 
more  pronounced  if  the  green  densities  had  been  ic^tical.  Ex- 
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Fie.  X  Scannine  electron  microeraphs  of  undopeU  alumina 
(97.5%  TD)  and  D-?  (99.0%  TD). 


tending  intermediate'Stage  sintering  leads  to  a  finer^gratned 
cenmie  in  the  case  of  alumina  containing  MgO  since  it  does 
not  undergo  the  same  degree  of  exaggerated  grain  and  pore 
growth  as  the  uiKloped  material.  The  high-resolution  SAXS 
experiments  were  performed  in  order  to  measure  the  Final- 
stage  samples  with  higher  precision  than  is  possible  with  the 
neutron  scattering  technique  and  to  extend  the  study  of  the 
evolution  of  the  pore  size  distribution  to  sample  densities 
greater  than  99%  TD. 

The  SAXS  Guinier  results  (Table  II  and  Fig.  4)  confirmed 
the  neutron  scattering  results  and  demonstrated  that  the  aver- 
agt  pore  radius  continues  to  increase  in  the  final  stage  of  sin- 
teri^  in  both  the  undoped  and  doped  alumina.  Furthermore, 
the  trend  continues  even  as  densities  reach  more  than  99'> 
TD.  When  these  results  are  combined  with  those  obtained 
from  the  FOcod  analysis  (Table  III  and  Fig.  5).  it  is  evident 
that  the  pore  size  distribution  broadens  with  sample  density 
in  both  eases. 

The  mote  qjuantitative  representation  of  the  pore  volume 
distributions,  obtained  by  the  maximum  entropy  technique 


13-7  (97.0%  TD).  13-8  (99.«q-  TD)  and  of  MgO-doped  samples  D-6 


(Fig.  6),  can  be  directly  related  to  the  grain  size  distributions 
observed  in  the  SEM  micrographs  (Fig.  7).  While  the  grain 
size  distribution  can  be  seen  directly  from  the  micrographs, 
the  SAXS  measurements  arc  used  to  obtain  a  statistically  sig¬ 
nificant  description  of  the  pore  size  distribution.  The  average 
pore  diameters  determined  from  the  pore  size  distribution 
(Table  IV)  quantitatively  confirm  the  neutron  seattering 
results  provided  that  the  volumetric  densities  are  used  in 
both  cases.' 

Earlier  studies'  in  which  the  dihedral  angles  between 
grains  were  measured  and  related  to  the  grain  size  distribu¬ 
tions  in  —99^  dense  alumina  and  alumina  with  0.25  wt% 
MgO  demonstrated  that  the  dihedral  angle  size  distribution, 
and  thus  the  grain  size  distribution,  were  narrower  for  the 
MgO-doped  sample.  However,  the  median  value  of  the  dihe¬ 
dral  angle  remained  the  same  for  both  undoped  and  MgO- 
doped  samples.  This  behavior  can  be  seen  in  the  micrographs 
for  the  ~97%  TD  samples  as  well  as  the  -99%  TD  sampto, 
although  the  micrographs  indicate  that  the  average  grain  size 
is  larger  for  the  undoped  cases.  The  narrowing  of  the  grain 
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Fit  8.  Results  from  earlier  multiple  small-angle  neutron  scatter¬ 
ing  studies'  showing  the  effective  pore  radius  as  a  function  of  sam¬ 
ple  density  for  the  undoped  (13-series)  and  MgO-doped  (D-series) 
samples. 

size  distribution  implies  that  the  ranges  of  the  surface  energy 
and/or  the  grain  boundary  energy  become  narrower  with  ad¬ 
dition  of  MgO. 

A  similar  trend  is  expected  for  the  pore  size  distribution  at 
the  same  densities.  The  derived  volume  distributions  for  the 
—99%  TD  samples  (13-K  and  D-7)  arc  plotted  together  in 
Fig.  9.  The  figure  suggests  that  the  size  distribution  is  indeed 
narrower  for  the  MgO-doped  sample  and  that  the  doped 
sample  contains  a  smaller  volume  of  pores  greater  than 
-SCO  nm  in  diameter.  Table  IV  supports  this  conclusion, 
showing  that  both  the  mean  pore  radius  and  its  standard  de¬ 
viation  are  smaller  for  the  sample  containing  MgO.  This  ef¬ 
fect  may  be  masked  somewhat  since  the  green  density  of  the 
MgO-doped  sample  was  higher  than  that  of  the  undoped 
sample,  meaning  that  the  average  pore  radius  of  the  MgO- 
doped  sample  was  larger  from  the  beginning.  If  both  samples 
had  started  ou*  'th  identical  green  densities,  the  differences 
observed  here  in  mean  pore  radius  and  in  the  width  of  the 
size  distributions  could  possibly  be  larger. 

The  pore  volume  fractions,  6.  calculated  from  the  maxi¬ 
mum  entropy  procedure  and  from  the  invariant,  where  the 


Fl^8.  Maximum  entropy-derived  volume  distributions  1^(0)  vs 
diameter,  D.  plots  for  ihe  -99%  dense  samples  (A)  13-8  (undoped) 
and  (B)  D-7  (MgO  doped). 


integral  in  Eq.  (4)  is  taken  from  <7  =  0  to<7m„  are  compared  to 
the  volumetric  determination  in  Table  V.  The  value  of  <l>  de¬ 
termined  from  maximum  entropy  is,  in  some  cases,  lower 
than  the  volumetric  determination,  which  has  an  uncertainty 
of  2:().003.  Furthermore,  the  value  determined  from  the  in¬ 
variant  generally  falls  above  that  determined  from  maximum 
entropy.  In  the  maximum  entropy  procedure,  an  analytical 
calculation  of  the  volume  fraction  based  on  the  known  scat¬ 
tering  function  for  spherical  scatterers  is  made.  On  the  other 
hand,  the  integral  used  in  calculating  the  invariant  begins 
at  17  0.  Since  no  data  are  available  at  values  below  q  = 

5  X  10"’  nm"',  {(q)  is  extrapolated  to<7  =  0  using  the  Guinier 
approximation  represented  by  Eq.  (1). 

Although  both  the  maximum  entropy  and  the  invariant  de¬ 
terminations  of  <t>  yield  volume  fractions  that  agree  well  with 
those  measured  volumetrically,  it  is  still  possible  that  a  small 
number  of  pores  in  each  sample  arc  too  large  to  be  measured, 
even  with  the  high-resolution  SAXS  instrument.  This  was  al¬ 
ready  suggested  from  the  Guinier  plots  (Fig.  2),  which  display 
noticeable  curvature  at  the  smallest  measured  q  values.  In  ad¬ 
dition,  the  Rf  values  obtained  from  the  Guinier  analyses 
(Table  11)  indicate  that  qR^  —  1.5  at  <7  =  5  x  lO"’  nm''. 
whereas  the  q  range  should  be  centered  upon  qR^  s.  1.2.'' 
Thus,  those  pores  with  sizes  that  correspond  to  the  extremes 
of  the  <7  range  may  be  underrepresented  in  the  pore  size  dis¬ 
tributions.  Further  evidence  of  the  existence  of  larger  pores  is 
provided  by  the  SEM  micrographs  in  which  a  few  pores  with 
diameters  larger  than  1  ^im  can  be  seen.  Except  that  the 
pores  at  the  extremes  of  the  q  range  may  not  be  completely 
included,  the  small  number  of  such  large  pores  will  not  signif¬ 
icantly  affect  the  derived  pore  size  distributions. 

VI.  Concluding  Remarks 

The  evolution  of  the  pore  size  distribution  during  final- 
stage  sintering  of  alumina  and  of  alumina  doped  with 
0.25  wi%  MgO  was  followed  using  high-resolution  small- 
angle  X-ray  scattering.  The  pore  radius  of  gyration,  obtained 
from  a  Guinier”  analysis,  was  observed  to  increase  from  -0.3 
to  —0.4  mhi,  independent  of  the  presence  of  MgO  sintering 
aid,  as  the  sample  density  increased  from  —91%  to  >99%  of 
theoretical  density  (Table  11).  In  contrast,  the  pore  radius 
derived  from  a  Porod''  analysis  and  assuming  a  spherical  pore 
shape  was  observed  to  decrease  from  -0.2  to  -0.1  >xm  over 
the  same  range  of  sample  densities  (Table  111).  Taken  together, 
these  results  indicate  qualitatively  that  the  pore  size  distribu¬ 
tion  broadens  with  increasing  sample  density.  In  addition, 
pore  volume  distributions,  which  provide  a  quantitative  de¬ 
scription  of  the  pore  size  distribution,  were  derived  for  each 
sample  using  a  maximum  entrooy  analysis.'"  Both  the  mean 
pore  radius  and  its  standard  deviation,  obtained  directly  from 
the  pore  size  distribution,  were  found  to  increase  with  in¬ 
creasing  sample  density.  The  trends  observed  here  confirm 
earlier  multiple  small-angle  neutron  scattering  studies'  * 
which  demonstrated  that  the  average  pore  radius  continues  to 
increase  in  the  final  stage  of  sintering  in  both  the  undoped 
and  the  doped  samples  and  suggested  that  this  trend  contin¬ 
ues  even  as  densities  reach  more  than  99%  of  theoretical. 
Furthermore,  the  average  pore  radii  determined  from  the 
SAXS  analyses  agreed  quantitatively  with  those  determined 
in  the  neutron  scattering  experiment,  provided  that  the  volu¬ 
metrically  determined  sample  densities  are  used  in  both 
cases. 

Scanning  electron  micrographs  of  both  undoped  (13-7  and 
13-S)  and  MgO-doped  (D-6  and  D-7)  samples  show  that  the 
majority  of  pores  are  smaller  than  I  /xm.  supporting  the  scat¬ 
tering  results.  The  average  grain  size  as  well  as  the  average 
pore  size  increases  for  both  doped  and  undoped  samf^es  as 
the  density  increases  from  97%  to  >99%  of  theoretical.  For  a 
given  sample  density,  the  average  grain  size  is  markedly 
smaller,  and  the  grain  size  distribution  is  narrower,  in  the 
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Table  V.  Scattering  Volume  Fractions 


Volume  fraction" 

Sample 

Volumetric 

measurement 

Maximum  entropy 
determination 

Invariant 

determination 

Undoped  alumina 

13-7 

0.03 

0.027 

0.029 

13-8 

0.01 

0.013 

0.012 

8-10 

0.023 

0.02 

0.025 

8-8 

0.014 

0.011 

0.015 

8-14 

0.007 

0.008 

0.010 

Doped  alumina 

D-6 

0.025 

0.023 

0.023 

D-7 

0.01 

0.011 

0.011 

D-7A 

0.003 

0.005 

0  007 

'The  unceriainty  in  all  pore  volume  fraction  values  is  ±0.0(J3. 


MgO-dopcd  samples.  This  confirms  earlier  studies^  which 
showed  that  the  distribution  of  dihedral  angles  between 
grains,  and  thus  the  grain  size  distribution,  in  "99%  dense 
MgO-doped  alumina  is  narrower  than  in  undoped  alumina  at 
the  same  density. 

Acknowledgment:  We  are  grateful  to  UKAEA  Harwell  for  the  MAXE 
software  used  to  perform  the  maximum  entropy  data  analysis. 
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Abstract 

Multiple  small-angle  neutron  scattering  was  used  to 
follow  the  evolution  of  the  pore-size  distribution  in 
a-Al20j  through  the  intermediate  and  final  stages  of 
sintering.  This  new  technique  enables  the  study  of 
microstructure  in  the  0-08-10  jim  size  regime,  which 
is  the  size  range  of  importance  for  many  materials 
systems,  without  needing  to  increase  the  resolution 
of  currently  available  small-angle  scattering  instru¬ 
ments.  The  microstructure  evolution  results  indicate 
a  nearly  constant  effective  pore  radius  for  the  alumina 
throughout  the  intermediate  sintering  stage,  ranging 
from  0-19  fimat  54%  of  theoretical  density  to  0-1 7  ^tm 
at  79%  dense.  As  the  alumina  densifies  further,  there 
is  a  transition  region  after  which  the  effective  pore 
radius  grows  rapidly  to  sO-fijim  at  97-5%  dense. 

1.  lotroductioo 

Small-angle  neutro.  'attering  (SANS)  is  an  impor¬ 
tant  technique  for  th  'tatistically  significant  evalu¬ 
ation  of  bulk  microstructural  parameters  in  materials 
containing  dilute  concentrations  [where  ‘dilute*  is 
represented  by  <p/(l  -<p)  —  <p  and  ^  is  the  scattering 
volume  fraction]  of scatterers  between  10  and  1000  A. 
However,  the  relevant  struaures  in  many  physical 
systems,  such  as  partially  sintered  ceramics,  can¬ 
not  readily  be  studied  by  SANS  since  the  usual 
approximations  are  not  valid  for  the  larger  sizes  and 
more-dense  structures  present  in  these  materials.  To 
perform  neutron  scattering  characterization  of  high 
technology  materials  as  a  function  of  thermal  process¬ 
ing,  it  has  become  increasingly  important  to  extend 
the  range  of  particle-size  determinations  usually 
associated  with  SANS  up  to  —lOjim  and  to  be  able 
to  include  the  analysis  of  thicker  and  less-dilute 
-  <p)>  <pj  systems  in  which  the  coherent  elastic 
differential  scattering  cross  section  is  dominated  by 
multiple  scattering.  These  issues  have  been  success¬ 
fully  addressed  by  the  multiple  small-angle  neutron¬ 
scattering  (MSANS)  formalism  of  Berk  &  Hardman- 
Rhyne  (1985,  1988),  which  applies  when  multiple 
scattering  arises  solely  from  the  compounding  of 
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single-particle  scattering  events,  meaning  that  there 
is  no  scattering  component  due  to  correlations 
between  particles.  Berk  &  Hardman-Rhyne  refer  to 
this  type  of  scattering  process  as  ‘incoherent  multiple 
scattering’. 

The  central  assumptions  of  the  MSANS  formalism 
are  that  the  observed  MSANS  intensity  is  the  result 
of  incoherent  multiple  scattering  from  a  randomly 
distributed  monodisperse  population  of  spheres  com¬ 
prised  of  material  of  uniform  density  (i.e.  containing 
no  internal  struaure)  and  that  the  phase  shift  u  that 
a  plane  wave  undergoes  in  traversing  a  particle  of 
radius  R  is  within  the  range  —0-1  to  2.  Effective 
(volume-weighted  average)  radii  of  the  scatterers  arc 
extracted  by  measuring  the  wavelength  dependence 
of  the  broadening  of  the  scattering  curves. 

An  earlier  study  of  microstruaure  evolution  in 
porous  silica  as  a  funaion  of  thermal  processing 
(Long  &  Krueger,  1989)  demonstrated  that  MSANS 
and  SANS  can  effectively  be  used  together  to  cover 
the  full  range  of  relevant  microstructure  sizes.  In  this 
work,  MSANS  was  used  to  follow  the  evolution  of 
the  pore-size  distribution  in  a-Al203  through  the 
intermediate  and  final  stages  of  sintering. 

Knowledge  of  the  microstructure  evolution  as  a 
function  of  thermal  processing  is  important  for  the 
development  of  process  models  in  ceramics.  The  pro¬ 
cess  models  for  both  the  silica  (Kingery  &  Berg,  1955) 
and  alumina  (Coble,  1961)  systems  are  comparatively 
well  understood.  Whereas  porous  silica  is  a  glassy 
system  which  sinters  by  means  of  viscous  flow, 
alumina  is  a  polycrystalline  system  which  sinters  by 
means  of  volume  and  surface  diffusion.  The  different 
sintering  mechanisms  should  give  rise  to  different 
microstructure  evolution.  Thus,  these  well  character¬ 
ized  materials  can  serve  as  model  systems  for  the  new 
MSANS  technique. 


2.  MSANS  theory 

The  interaction  of  neutrons  with  matter  is  character¬ 
ized  by  the  phase  shift  u  that  a  plane  wave  undergoes 
in  traversing  a  particle  of  radius  R.  This  phase  shift 
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determines  the  shape  of  the  single-particle  differential 
scattering  cross  section,  d2(Q)/df2,  as  a  function  of 
the  scattering  wavevector  Q,  where  |Ql  =  ItteIK  and 
e  <  1  is  the  scattering  angle,  v  depends  upon  Ap,  the 
contrast  of  the  particle  or  void,  relative  to  the  scatter¬ 
ing  matrix  such  that 

v  =  2ApRX,  (1) 

where  Ap  is  written  as 

Ap=l(b,/V)-p^,,^,.  (2) 

celt 

The  bi  are  the  coherent  scattering  lengths  of  the 
individual  nuclei  whereas  p„.(rix  is  the  average  scatter¬ 
ing-length  density  of  the  matrix.  Thus,  the  phase  shift 
is  directly  dependent  upon  the  material  contrast  (Ap), 
the  dimension  of  the  scattering  particle  or  void  {2R) 
and  the  neutron  wavelength  (A). 

If  vxi  (the  usual  SANS  diftaction  regime),  the 
measured  scattered  intensity  can  be  described  by  the 
Bom  approximation  (Porod,  1951;  Weiss,  1951)  of 
scattering  theory  or  the  Rayleigh-Gans  model  (Ray¬ 
leigh,  1910,  1914;  Gans,  1925)  of  wave  optics.  For 
appropriately  thin  samples,  the  majority  of  the 
incident  neutron  beam  is  transmitted  through  the 
sample  and  the  remaining  neutrons  are  scattered  only 
once  before  exiting  the  sample.  The  scattered  intensity 
as  a  function  of  Q  is  independent  of  neutron 
wavelength  and  its  shape  near  Q  =  0  depends  only 


1(0) 


(a) 


1(0) 


Rg.  t.  (a)  ReprcMnuUon  of  an  I{Q}osQ  fcanered-intensity 
curve  mrising  from  convenlional  smaJI-angie  scattering.  iUiutrai- 
ing  its  wavelength*  independent  shape  and  the  transmitted-beam 
region  near  Q"0.  (h)  Representation  of  a  typical  set  of 
/(Q)  o$Q  multiple  small-angle  scattered-intensity  curves  from 
a  single  sample,  illustrating  their  wavelength-dependeni  shape 
and  the  empty-beam  curve. 


on  the  particle  dimensions.  Fig.  l(n)  represents  a 
typical  SANS  scattering  curve.  The  scattered  intensity 
around  Q  =  0  cannot  be  measured  directly  because  it 
lies  in  the  same  region  as  the  transmitted  beam,  which 
is  10^  to  10’  times  more  intense  than  the  scattered 
beam.  Therefore,  a  beamstop  is  usually  employed  to 
prevent  the  transmitted  beam  from  reaching  and  thus 
damaging  the  neutron  detector.  For  any  particle 
shape,  the  particle  size  can  be  described  by  its  radius 
of  gyration  R,  or  Guinier  radius  (Kostorz,  1979) 
which  is  the  same  as  the  radius  of  gyration  in 
mechanics  and  is  obtained  from  an  analysis  of  the 
small-Q  portion  of  the  scattered-intensity  curve.  The 
range  for  which  the  Guinier  analysis  is  valid  is  defined 
by  "“1-2  (Kostorz,  1979).  In  addition,  the  ratio 
of  the  total  scattering  surface  area  to  the  volume  of 
the  sample  in  the  neutron  beam,  S/V^,  can  be 
obtained  from  an  analysis  of  the  large-Q  portion  of 
the  scattered-intensity  curve.  In  this  region,  the  scat¬ 
tered  intensity  obeys  the  <?"*  Porod  law  (Kostorz, 
1979).  An  equivalent  spherical  radius  can  be  deter¬ 
mined  from  the  total  scattering  surface  area. 

At  the  opposite  extreme,  when  v  >  1  (the  multiple 
refraction  regime),  every  particle  scatters  and  the 
scattered  intensity  is  derived  from  geometrical  optics 
(von  Nardroff,  1926).  The  shape  of  the  scattering 
curve  is  independent  of  the  particle  radius  and 
depends  only  upon  the  difference  in  index  of  refrac¬ 
tion  between  the  scatterers  and  the  matrix  which,  in 
turn,  depends  upon  neutron  wavelength  as  A^  lliis 
regime  is  not  pursued  further  here. 

When  v^'l  (the  MSANS  diffraction  regime),  the 
scattered  intensity  cannot  be  described  by  either  con¬ 
ventional  SANS  diffraction  or  by  multiple  refraction. 
Scatterers  with  radii  in  the  0*  1-10  um  range  generally 
undergo  phase  shifts  of  the  order  of  1  for  the  range 
of  neutron  wavelengths  usually  available  for  SANS. 
Since  larger  particles  scatter  more  intensely  than 
smaller  ones,  the  likelihood  of  multiple  scattering 
becomes  greater  (Berk  &  Hardman-^yne,  19gg). 
Furthermore,  die  incidence  of  coherent  scattering 
between  larger  particles  is  greatly  decreased,  meaning 
that  the  single-panicle  scattering  cross  section  does 
not  contain  a  component  due  to  interpanicle  interfer¬ 
ence.  Thus,  in  the  MSANS  diffraction  regime,  multi¬ 
ple  scattering  arises  solely  from  the  statistical  com¬ 
pounding  of  numerous  coherent  elastic  single-particle 
scattering  events  in  which  interpanicle  interference 
does  not  play  a  significant  role  (Berk  &  Hardman- 
Rhyne,  1985). 

A  measure  of  the  amount  of  incoherent  multiple 
scattering  is  given  by  the  parameter  f  (Berk  &  Hard- 
man-Rhyne,  1988),  which  is  the  sample  thickness,  z, 
divided  by  the  statistical  neutron  mean  free  path 
length,  L 

z  =  z/I=l-5ipz(ApX)^R,  (3) 

where  ip  is  the  volume  fraction  of  scatterers.  z  is 
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related  to  the  sample  transmission  as  T,  =  exp  {—2}. 
Since  Ap  involves  only  coherent  scattering  lengths,  2 
can  only  be  calculated  from  the  measured  sample 
transmission  after  correcting  for  absorption  and  inco¬ 
herent  single- particle  scattering,  where  applicable.  A 
value  of  2  =  0-0  means  that  there  is  no  multiple  scatter¬ 
ing.  In  practice,  the  data  from  SANS  samples  of  finite 
thickness  such  that  0<2^]  always  contain  some 
multiple  scattering  which  is  either  negligibly  small  or 
treated  as  a  small  correction  to  the  single-particle 
scattering.  In  MSANS  measurements,  the  amount  of 
multiple  scattering  is  copious  and  2  typically  ranges 
from  10  to  400. 

In  the  region  of  validity  of  MSANS,  the  intensity 
of  the  transmitted  beam  is  immeasurably  small  and 
the  width  of  the  scattered  intensity  curve  near  Q  =  0 
is  broadened  far  beyond  the  broadening  due  to  instru¬ 
mental  resolution.  The  amount  of  beam  broadening 
is  dependent  upon  the  incident  neutron  wavelength 
as  illustrated  in  Fig.  1(h)  where  representative 
MSANS  curves  are  shown  for  different  wavelengths 
incident  on  a  single  sample.  The  intensities  have  been 
normalized  such  that  7(0)  =  1-0.  In  each  case,  the 
scattered  intensities,  7(Q),  have  a  curvature,  r„  near 
Q  =  0  and  can  be  approximately  described  as  a 
Gaussian  with  a  width  proportional  to  A^.  For  large 
Q,  there  is  a  crossover  to  single-scattering  Porod 
behavior  (Berk  &  Hardman-Rhyne,  1985),  where  the 
cross  section  is  proportional  to  Q~*. 

It  is  evident  from  (I )  and  (3)  that  the  contribution 
of  incoherent  multiple  scattering  to  the  total  cross 
section  depends  upon  a  balance  of  the  parameters 
which  define  v  and  f.  For  a  given  neutron-scattering 
contrast  and  wavelength,  (3)  indicates  that  the 
amount  of  multiple  scattering  depends  linearly  upon 
the  sample  thickness  and  the  radii  of  the  individual 
scatterers.  Therefore,  even  though  in  practice  thick 
samples  (0’S-l*0cm)  are  used  to  ensure  copious 
multiple  scattering,  the  MSANS  formalism  can  also 
be  applied  to  much  thinner  samples  (0-1 -0-2  cm)  if 
the  radii  of  the  scatterers  are  sufficiently  large.  (From 
a  practical  standpoint,  sample  thickness  can  be 
chosen  so  that  the  wavelength  dependence  of  the 
width  of  the  scattering  curve  is  easily  observed  within 
the  range  of  available  neutron  wavelengths.) 

Although  the  wavelength  dependence  at  small  Q 
seems  to  be  the  same  as  that  for  the  multiple  refraction 
regime,  the  underlying  mechanism  is  different.  When 
v^l,  an  effective  radius,  /ietj(0),  where  the  0  refers 
to  <?  =  0,  can  be  determined  for  the  scatterers  from 
the  wavelength  dependence  of  using  the  MSANS 
formalism  represented  in  equations  (2.12)-(2.IS)  of 
Berk  &  Hardman-Rhyne  (1985,  1988).  Tbe 
dimensionless  parameter  rc[<?7?3  =  where 

[ . . .  ]  means  with  respect  to,  is  related  to  2  via  an 
empirical  relation  (Berk  &  Hardman-Rhyne,  1988) 

rJOR]  =  0-926{f  In®  “  f}®  ’  for  2  >  5,  (4) 


which  approximates  the  theory  to  better  than  0-1% 
when  2  >  1 0. 

The  MSANS  formalism  in  effect  replaces  the  stan¬ 
dard  Guinier  analysis  in  the  the  small-Q  region  of 
/(Q)  where  R,ff(0),  rather  than  /?,,  defines  the  size 
of  the  scatterers.  When  polydispersity  is  present, 
can  be  written  in  terms  of  the  moments  of  an 
assumed  particle-size  distribution.  For  the  small-(? 
regime,  ( Berk  &  Hardman-Rhyne, 

1985),  where  the  brackets  denote  the  arithmetic 
average  of  the  argument.  Note  that  in  the  low- 2  limit 
the  MSANS  formalism  no  longer  applies  and  the 
moments  are  expressed  in  terms  of  the  Guinier  radius 
of  gyration  such  that  Rg  =  [3(R*>/5(i?®>]'*'‘  (Kostorz, 
1979).  If  a  similar  relation  for  the  effective  radius  can 
be  found  using  the  large- Q  regime,  then  the  average 
radius  as  well  as  the  width  of  the  size  distribution 
can  be  found  as  shown  below. 

The  standard  Porod  analysis  (Kostorz,  1979)  can 
be  applied  to  the  large-Q  region  of  HQ).  In  this 
region,  the  scattered  intensity  can  be  written  as 

I{Q)^PQ-*+B,  (5) 

where  P  is  Porod's  constant  and  B  is  a  background 
term.  The  normalized  total  surface  scattering  area  per 
unit  volume,  S/  VJ,,  is  determined  directly  from  P.  An 
effective  pore  radius,  ReB(QO),  where  00  refers  to  the 
large-<?  regime,  is  calculated  directly  from  S/ 
assuming  spherical  scatterers.  For  polydisperse  sys¬ 
tems,  the  moments  of  the  size  distribution  are 
expressed  as  Rt„(oo)-{R^)/{R^).  If  both  R^niO)  and 
Ren(<»)  are  known,  (R)  as  well  as  the  width,  P,  of  an 
assumed  particle-size  distribution  can  be  determined. 

To  apply  the  MSANS  formalism  to  the  analysis  of 
scattering  from  a  system  of  spheres  containing  no 
internal  structure,  (3)  and  (4)  indicate  that  the  sample 
thickness,  the  scattering  volume  fraction  (which  is  a 
refinable  parameter),  the  scattering-length  density 
and  the  wavelength  dependence  of  must  be  known 
in  order  to  derive  Reft(O).  In  addition,  the  constraints 
that  2^10  and  0-lsi/s2  must  be  satisfied.  When 
j'>2,  the  scattering  intensity  falls  into  a  crossover 
regime  between  diffraction  and  refraction.  However, 
if  2  is  sufficiently  large,  the  MSANS  formalism 
may  even  apply  if  i/>2  (Berk  &  Hardman-Rhyne, 
1988). 

The  assumption  that  the  scattering  system  is  domi¬ 
nated  by  incoherent  multiple  scattering  is  extremely 
important  since  it  means  that  the  scattered  intensity 
does  not  contain  an  interparticle  interference  term 
due  to  pair-pair  correlations  between  particles.  Since 
the  scatterers  are  so  large,  Le.  comparable  in  size  to 
the  mean  free  path  length  of  the  neutrons,  there  is 
no  coherent  scattering  between  particles.  Thus,  even 
a  dense  system  of  scatterers  may  be  considered  dilute 
from  the  perspective  of  MSANS  even  though  ip /(I  - 
ip) —  ip  does  not  apply. 
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Table  1.  Alumina  samples  prepared  from  Baikowski 
CR6,  6  g“'  powder 


Thickness 

Pore  volume 

%TD 

Sample 

(cm) 

fraction,  if 

Volumetric 

MSANS 

13-2 

1-019 

n*6 

.‘i6-5 

54 

13-3 

0-9SS 

0-35 

67-0 

65 

13-« 

0  g94 

0-21 

79-1 

79 

13-5 

0-914 

015 

85-4 

85 

13-1 

0-SSl 

O-II 

90-8 

89 

13-< 

0-852 

0-075 

94-6 

92-5 

13-7 

0-840 

0-055 

97-1 

94-5 

1 3-8 

0  804 

O-p* 

98-9 

97-5 

3.  Alumina  samples 

The  samples  were  prepared  from  Baikowski  alumina 
powder  (CR6,  6  m^  g"'),  slipcast  lo  an  average  green 
density  of  53%  of  theoretical  density  (TD),  where 
TD  =*  3-965  gcm“\  The  slip  consisted  of  43-5vol.% 
alumina  powder  mixed  into  water  along  with  0-5  wt% 
of  a  polyelectrolyte  dispersant  (Darvan  C).  The  sus¬ 
pension  was  ultrasonicated  for  15  min  and  allowed 
to  set  overnight  before  being  poured  into  molds  with 
Teflon  rings  on  top  of  plaster-of-paris  blocks.  After 
55  min,  the  molds  were  removed  from  the  blocks  and 
sealed  in  a  container  overnight.  The  partially  dried 
samples  were  then  removed  from  the  molds  and 
sealed  in  containen  for  2  d.  The  containers  were  then 
opened  to  the  air  for  an  additional  5-7  d  before  a 
series  of  oven-drying  steps  was  performed  as  follows: 
2  d  at  337  K,  1  d  at  354  K,  I  d  at  375  K  and  1  d  at 
389  K,  after  which  a  dry  solid  green  body  is  obtained. 
In  this  fashion  a  family  of  eight  nearly  identical 
samples  was  prepared. 

During  sintering,  samples  were  removed  one  by 
one  at  different  times  from  the  furnace  to  yield 
examples  of  material  between  56-5  and  97-5%  TD. 
The  charaaeristics  of  the  alumina  samples  are  listed 
in  Table  I.  Sample  densities  measured  volumetrically 
as  well  as  those  determined  from  the  MSANS 
measurements  are  given.  The  sintering  time- 
temperature  curve  is  shown  in  Fig.  2.  Each  circle 
represents  a  sample  removed  from  the  furnace  and 
the  number  alongside  each  circle  is  the  measured 
%TD. 

4.  Neutron'Scattering  measurements 

Measurements  were  made  at  the  20  MW  research 
reactor  at  the  National  Institute  of  Standards  and 
Technology  (NIST).  The  NIST  SANS  facility 
(Glinka,  1981;  Glinka,  Rowe  &  LaRock,  1986)  makes 
use  of  a  velocity  selector  to  choose  the  mean  incident 
wavelength  of  the  neutrons.  A,  which  have  a 
wavelength  spread  ^A/A=0'25.  Copious  neutron 
fluxes  can  be  obtained  in  the  wavelength  region 
between  5  and  20  A  due  to  a  cryogenic  moderator 
installed  in  the  reactor  core.  The  neutron  beam  is 
collimated  with  a  12  mm  aperture  which  follows  the 
velocity  selector  and  an  8  mm  aperture  which  pre¬ 


cedes  the  sample  position.  Scattered  neutrons  are 
detected  by  a  64  x  64  cm  position-sensitive  detector 
divided  into  128  x  128  pixels.  The  data  are  circularly 
averaged  to  produce  one-dimensional  intensity,  HQ), 
vs  Q  curves.  Intensities  measured  for  Porod  analysis 
were  placed  on  an  absolute  scale  by  comparison  with 
HQ)  from  a  calibrated  silica-gel  sample  measured 
under  the  same  conditions. 

The  MSANS  measurements  were  performed  with 
no  beamstop  in  front  of  the  detector  since  the  intense 
transmitted  beam  associated  with  conventional  SANS 
measurements  is  absent.  The  absence  of  the  beamstop 
ensures  that  the  beam  center  can  be  found  with  high 
precision.  With  the  center  of  the  Gaussian  precisely 
known,  an  accurate  determination  of  the  radius  of 
curvature  can  be  made.  Measurement  times  ranged 
from  10  min  to  1  h,  depending  upon  the  neutron 
wavelength.  At  the  shortest  wavelengths,  the  incident 
neutron  beam  is  attenuated  as  needed  in  order  to 
avoid  saturating  the  detector.  For  a  given  wavelength, 
the  broadest  scattered  intensity  curves  require  the 
longest  measuring  times  in  order  to  accumulate 
sufficient  counts. 

5.  Data  analysis 
5.1.  MSANS  analysis 

MSANS  measurements  were  performed  on  each 
sample  at  a  minimum  of  five  neutron  wavelengths. 
The  scattering  curves  exhibited  the  expected  broaden¬ 
ing  with  increasing  neutron  wavelength  as  shown  in 


Fig.  2.  Sintering  time-iempeniure  curve  for  the  alumina  samples. 
Each  circle  represents  a  sample  removed  from  the  Turaace  and 
the  number  alongside  is  the  percent  of  theoretical  density. 
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Fig.  3,  where  1{Q)  vs  Q  curves,  normalized  to  /(O)  = 
1-0,  are  plotted  for  the  85%  TD  sample.  The  features 
at  Q^sO'OOSA'',  most  obvious  for  the  longcr- 
wavelength  data,  are  due  to  a  change  in  deteaor 
sensitivity  at  the  beam-center  position.  To  avoid  this 
region,  Gaussian  fits  were  made  to  each  curve  in  the 
Q  range  beginning  where  the  intensity  is  —85-90% 
of  /(O)  and  extending  to  —45%  of  /(O).  Below  40% 
of  1(0),  the  scattering  curves  did  not  fit  well  to  a 
Gaussian. 

For  each  curve  fit,  a  radius  of  curvature  is  extracted 
from  the  standard  deviation,  a,  such  that  rc  = 
Blank  spectra,  measured  at  each 
wavelength  with  no  scatterer  in  place,  are  also  fitted 
to  enable  a  correction  for  the  instrument  function 
(Hardman- Rhyne  &  Berk,  1985).  Radii  are  extracted 
in  each  case  from  the  empirical  formula  in  (4)  using 

(3)  for  z  These  values  are  averaged  to  obtain  an 
effective  radius,  i?cfi(0),  for  the  scatterers  in  each 
sample. 

Finally,  both  i?,(i(0)  and  the  scattering-volume 
fraction,  <p,  are  refined  simultaneously  for  each 
sample  using  a  least-squares  fitting  method.  Owing 
to  the  approximate  nature  of  the  empirical  formula 

(4) ,  the  full  numerical  MSANS  analysis  describing 
as  a  function  of  A,  developed  in  (2.I2)-(2.I5)  of  Berk 
&  Hardman- Rhyne  (1985),  is  used  in  the  refinement 
procedure. 

5.2.  Porod  analysis 

The  scattering  at  large  Q  was  measured  for  each 
sample  with  an  incident  neutron  wavelength  of  6  A 
and  the  detector  offset  to  an  angle  of  3*5°.  Sample 
transmissions  were  determined  by  dividing  the  total 
detector  counts  accumulated  for  each  sample  by  the 


R|.  3.  Measured  MSANS  /(Q)  vjQ  curves  from  the  85%  TD 
alumina  sample  for  ei|ht  neutron  wavelengths.  Die  curve  Tor 
the  empiy-beam  ‘blank*  was  taken  at  A  >  13  A. 


total  deteaor  counts  accumulated  under  the  same 
conditions  with  no  sample  in  place.  This  method 
yields  the  equivalent  of  a  single-scattering  sample 
transmission.  Values  between  0-87  and  0-97  were 
measured,  which  fall  within  the  range  expeaed  for 
AI2O3  when  no  multiple  scattering  is  present. 

Since  the  scattered  intensity  follows  (5)  in  the 
Porod  regime,  both  the  Porod  constant,  P,  and  the 
background  term,  B,  were  found  direaly  by  making 
linear  least-squares  fits  to  IQ*  vs  Q*  curves  for  each 
sample.  The  largest  possible  linear  region  of  the 
curves  was  fit.  In  each  case,  the  maximum  Q  value 
was  0-16  A"',  which  is  the  largest  Q  value  for  which 
data  were  obtained.  The  minimum  Q  value  ranged 
from  0-035  A“'  for  the  97-5%  TD  sample  to  0-1  A~'  for 
the  54%  TD  sample.  The  ratio  of  the  normalized 
total  scattering  surface  area  to  the  volume  of  sample 
in  the  neutron  beam  S/  V*  is  determined  direaly  from 
P.  With  the  assumption  of  spherical  scatterers, 
for  the  large-G  regime  is  also  determined. 


6.  Results  and  discussion 

The  measured  values  of  (symbols)  are  plotted 
versus  A  in  Fig.  4  along  w’th  the  theoretical  values 
(solid  lines)  determined  using  the  MSANS  formalism 
[(2.I2)-(2.15)  in  Berk  &  Hardman-Rhyne  (1985)].  As 
the  scattering  volume  fraaion  decreases  {Le.  the  %TD 
increases),  the  calculated  values  of  K,fl(0)  become 


Fig.  4.  Measured  radii  or curvaiure  (symbols)  and  the  correspond¬ 
ing  MSANS  theoretical  fit  (solid  lines)  as  a  function  of  neutron 
wavelength.  The  sample  densities  are  those  refined  from  the 
neulron-scattering  measurements. 
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Table  2.  MSANS  results 


Neutron 
wavelength 
A  (Al 


Sample 

%TD 

Min. 

Max. 

13-2 

34 

6  0 

9  9 

13-3 

65 

6  0 

9-9 

13-4 

79 

6  0 

10  9 

13-5 

85 

6  0 

130 

I3.| 

89 

7  0 

14  2 

13-6 

92  5 

9  9 

14  2 

13-7 

94  5 

9  9 

16  5 

13-8 

97  5 

10  9 

16  5 

Pore  radius 
R„(0)(»im) 

Min. 

Max 

0-20<l) 

152 

414 

0  19(1) 

98 

289 

0- 170(5) 

49 

193 

O' 1811) 

32 

204 

O' 23(2) 

47 

234 

0-32(3) 

96 

298 

0-53(5) 

73 

248 

0-75(15) 

40 

126 

Ratio 

■ 

Phase  shift 

y 

Min. 

Max. 

Min 

Max 

14 

26 

0  13 

0  22 

1 1 

21 

0-12 

0  21 

7 

17 

0  10 

0  22 

6 

17 

0  09 

0-28 

7 

19 

0  (7 

0-41 

II 

21 

0  35 

0  66 

9 

19 

0  45 

0  93 

6 

13 

0-33 

1  10 

*  Ratio  =  it,«(0)/ R,(0).  where  is  the  pore  radius  calculated  assuming  no  muhiple  scattering  is  present. 


much  more  sensitive  to  small  changes  in  the  volume 
fraction,  <p.  For  example,  an  uncertainty  of  0-01  on 
the  value  <p  =  O-OI  has  a  much  greater  effect  on  /{efl(O) 
than  the  same  uncertainty  on  the  value  <p  =  0'4.  In 
these  cases,  the  refinement  procedure  was  especially 
useful  as  it  allowed  the  simultaneous  determination 
of  (p  and  Rt(j(0)  values  which  best  fit  the  r,  vs  A  data. 
Values  for  the  phase  shift  v  ranged  from  0-1  to  LO, 
as  shown  in  Table  2,  which  fall  well  within  the  allowed 
range  for  MSANS  analysis. 

The  slopes  of  the  w  A  curves  are  steepest  for  the 
intermediate-stage  alumina  samples,  which  have  the 
largest  pore  volume  fraction.  These  samples  also  have 
the  largest  z  values  and  thus  the  largest  amount  of 
multiple  scattering,  even  though  the  effective  pore 
radii  are  smaller  than  those  of  the  final-stage  samples. 
In  addition,  the  scattering  intensities  show  beam 
broadening  at  smaller  A  values  and  the  width  of  the 
Gaussian  at  Q  =  0  increases  more  rapidly  with 
increasing  A. 

When  the  pore  volume  fraction  is  small,  as  in  the 
97-5%  TD  sample,  f  is  correspondingly  smaller  and 
the  slope  of  the  r^vsk  curves  is  much  flatter.  Since 
the  amount  of  multiple  scattering  is  much  less  in  the 
denser  samples,  the  scattered  intensities  must  be 
measured  using  longer-wavelen'gth  neutrons  for  beam 
broadening  to  be  observed.  Without  the  cryogenic 
moderator,  copious  neutron  fluxes  at  wavelengths 
beyond  —12  A  would  not  be  available  and  the  most 
dense  samples  could  not  have  been  measured  using 
the  MSANS  technique. 

The  refined  value  for  the  effective  pore  radius 
obtained  from  the  data  in  Fig.  4,  is  listed  for 
each  sample  in  Table  2,  along  with  the  ranges  of  A, 
i  and  v.  Also  listed  is  the  ratio  of  Rn^O)  to  the 
equivalent  spherical  radius  i{,(0)  calculated  assuming 
that  no  multiple  scattering  is  present,  i.e.  using  the 
standard  Guinier  analysis  to  obtain  and  then 
assuming  spherical  scatterers  to  obtain  R,(0). 

The  data  points  from  Fig.  4  for  all  of  the  samples 
(circles)  along  with  the  empirical  curve  as  given  by 
(4)  (solid  line)  are  plotted  as  a  function  of  f  in  Fig. 
S.  The  dimensionless  parameter  rc[Qi{l  is  obtained 
by  multiplying  r,[Q]  by  the  best  R^O)  value  for  each 


sample  after  refinement.  The  f  values  for  these 
samples  range  from  40  to  400,  well  within  the  range 
of  validity  of  the  empirical  formula.  Thus,  excellent 
agreement  with  (4)  is  expected. 

A  In  (/)  vs  In  (0)  plot  of  the  Porod  data  (circles) 
for  the  85%  TD  sample  is  shown  in  Fig.  6.  The  fitted 
background  term  [see  (5)]  of  0-0l8(  1)  cm'*  was  sub¬ 
tracted  prior  to  plotting.  The  solid  line  is  the  linear 
fit  to  ln(/)  vs  ln(0)  in  the  range  0-055s(?s 
0-155A”‘,  where  the  slope  was  determined  to  be 
-4-02(2).  The  effeaive  pore  radius,  and  the 

minimum  Q  value  used  for  the  Porod  fit  are  listed  in 
Table  3  for  each  sample,  along  with  the  normalized 
total  scattering  surface  area,  5/  Vj,.  The  normalized 
pore  volume  fraaion  ip/iS/V,,)  and  the  number  of 
pores  per  unit  sample  volume,  N/  Vj,,  are  also  tabu¬ 
lated  for  each  sample. 

The  Re«(0)  and  R*n(<»)  values  calculated  from  the 
MSANS  and  Porod  measurements,  respectively,  are 
plotted  in  Fig.  7,  along  with  those  obtained  using 


Fi|.  S.  Combined  plot  of  rJQft]  et  I  for  all  saniplet  ai  all  neutron 
wavelengths  measured.  Tlic  solid  line  u  the  empirical  formula 
(4).  The  data  points  are  the  products  of  the  measured  r,{Q]  and 
the  best-fit  value  for  the  eflective  pore  radius,  for  each 

sample. 
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mercury  intrusion  porosimetry  (Rootare  &  Prenzlow, 
1967)  for  the  samples  below  85%  TD,  where  open 
porosity  is  expected.  The  MSANS  determinations, 
which  are  weighted  to  the  larger  volume  scatterers 
(Long  &  Krueger,  1989),  indicate  that  the  average 
pore  radius  for  the  distribution  within  the  alumina  is 
essentially  flat,  with  a  value  of  '-0-2  ^jim,  in  the  inter¬ 
mediate  sintering  stage  and  increases  sharply  to 
aO-6  p.m  in  the  final  stage.  The  large  error  bar  on  the 
value  for  the  97-5%  TD  sample  is  a  consequence  of 
the  shallow  slope  in  the  vs  A  curve  (Fig.  4).  Since 
three  local  minima  in  the  fitting  parameter  were 
found  for  values  between  0-6  and  0-9  p-m,  no 

single  fit  could  be  chosen  as  the  ‘best  fit'  to  the 
MSANS  formalism.  An  uncertainty  in  the  pore 
volume  fraction,  <p,  of  approximately  O-OOS  was  also 
found  from  the  fitting  procedure. 

The  /?,B(ac)  values  determined  from  the  Porod 
analysis  are  weighted  to  the  particles  with  the  smallest 
scattering  surface  areas.  Although  consistently  lower, 
the  values  agree  qualitatively  with  the  RtgiO) 

values.  The  average  pore  radius  remains  nearly  con¬ 
stant  at  a  value  of  —0- 1  p.m  throughout  the  intermedi¬ 
ate  stage  and  increases  to  aO-4  p.m  in  the  final  stage. 
The  large  error  bar  on  Rta{°o)  for  the  97-5%  TD 
sample  is  mainly  due  to  the  uncertainty  in  «p  as 
determined  from  the  MSANS  analysis. 

The  porosimetry  results  agree  qualitatively  with 
results  obtained  from  both  the  MSANS  and  Porod 
analyses  for  intermediate-stage  sintering.  The 
observed  quantitative  discrepancies  are  expected 
since  the  MSANS  and  Porod  radii  are  derived  from 
a  measure  of  all  of  the  voids  in  the  system,  with 
MSANS  radii  weighted  to  the  larger-volume  and 
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Rf.  6.  Pr^rod  dau  Tor  ihe  85%  TD  sample  ploued  on  a 
In  (/)  M  In  (Q)  scale  after  subtraction  of  the  6ued  background 
term  [see  (3)].  The  solid  line  is  the  linear  least-squares  fit  to  the 
dau  in  the  range  0-055S  (?s0'l55  A"'.  The  fitted  slope  is 
-4-02(2). 


Porod  radii  weighted  to  the  smaller-volume  voids.  On 
the  other  hand,  the  porosimetry  radii  mainly  represent 
the  size  of  the  small  necks  between  pores. 

Mercury  intrusion  porosimetry  assumes  that  the 
largest  po;  “s  are  on  the  outside  of  the  sample.  This 
is  important  because  the  pressure  required  to  fill  a 
large  pore  is  less  than  that  required  to  fill  a  small 
pore.  However,  as  the  mercury  peneirates  a 
homogeneous  sample  consisting  of  a  random  distri¬ 
bution  of  necks  and  larger  channel-like  voids,  a  sig¬ 
nificant  number  nf  necks  may  be  encountered  before 
the  mercury  reaches  the  open  channels.  Since  the 
neck  size  is  measured  at  a  higher  pressure  than  that 
required  for  the  open  channel,  the  mercury  reaches 
the  channel  after  the  appropriate  pressure  to  measure 
its  size  has  been  exceeded.  Thus,  large  voids  in  the 
interior  of  a  sample  may  not  be  measured,  even  if 
they  form  a  majority  population. 

The  MSANS  and  Porod  results  for  pore  size  a&iume 
spherical  scatterers  whereas  the  porosimetry  results 
are  obtained  assuming  cylindrical  pore  morphology. 
Yet  the  same  trend  in  the  evolution  of  the  effective 
pore  radius  during  intermediate  stage  sintering  is 
reported  in  each  case.  This  occurs  because  the  average 
length  of  the  cylindrical  channels  is  ^10  pm,  which 
is  beyond  the  measurable  range  of  even  the  MSANS 
technique.  Thus,  the  pore  channels  can  be  effeaively 
represented  by  strings  of  spheres,  with  RtniO)  and 
corresponding  to  the  cross-sectional  radii  of 
the  cylinders.  /?efl(0)  is  dominated  by  the  larger- 
diameter  and  Reii(oo)  by  the  smaller-diameter  chan¬ 
nels  in  the  size  distribution. 

The  important  result  of  the  MSANS,  Porod  and 
porosimetry  measurements  is  that  the  effective  pore 
radius  remains  constant  through  the  intermediate 
stage  of  sintering.  This  is  consistent  with  topological 
decay  models  of  the  intermediate  sintering  stage  in 
which  an  interconnected  pore  network  decays  in  a 
stable  manner  (Rhines  &  DeHoff,  1984).  Throughout 
this  stage,  the  network  of  interconnected  channels 
decays  such  that  the  ratio  v>/{S/  V^),  or  pore  volume 
fraaion  to  surface  area  per  unit  volume  of  sample  in 
the  neutron  beam,  remains  constant  (DeHoff,  Rum- 
mel,  LaBuff  &  Rhines,  1966).  Consequently,  the  pore 
channels  become  fewer  while  retaining  a  constant 
diameter  as  densification  proceeds.  This  can  easily 
be  observed  in  Tables  2  and  3,  which  show  constant 
R»fl(0),  R.fl(co)  and  <p/{SfvJ)  values  for  the  inter- 
mediate-sintering-stage  samples.  .-Mso  shown  is  the 
number  of  pores  per  unit  sample  volume  (N/V^), 
which  decreases  sharply,  as  expected,  as  intermedi¬ 
ate-stage  sintering  progresses. 

Before  the  channels  have  disappeared  altogether, 
a  transition  into  the  final  sintering  stage  occurs  in 
which  the  porosity  becomes  isolated.  Pore  isolation 
occurs  because  the  channels  rupture  or  ‘pinch  off 
when  they  can  no  longer  maintain  their  length  while 
at  the  same  time  remaining  stable  against  the  increas- 

38 


S  KRUEGER.  G.  G.  LONG  AND  R.  A.  PACE 


289 


Table  3. 

Porod  results 

Normalized 

.Number  of 

Normalize': 

pore  -. olume 

pores  per 

surface  area 

fraciion 

unit  '.olume 

Pore  radius 

s\;’ 

ai/fS/K) 

V.  K 

Sample 

%TD 

')  lUi,(x)  (tim» 

1  xi0*cni"‘i 

(xiO  ‘cm) 

( X lO'"  cm"'*i 

13-2 

54 

0  10 

0  09i(i> 

if'2(2 

3-0i(ai 

)4'?  01 27) 

13-3 

65 

0  09 

0  1O811 1 

9-7(1) 

3  60(4) 

(16-0(5) 

13-4 

79 

0  07 

0  108I11 

5-83(3) 

3  60(2) 

39  6(48) 

13-5 

85 

0  065 

0  120(  1  > 

3-74(2) 

4  01(2) 

20  5(18) 

13-1 

89 

0  055 

0  154(11 

2-15(1) 

5)2(2) 

7  2(10) 

13-6 

92  5 

0  045 

0-185(11 

1-22(1) 

6- 16(3i 

2  84(72) 

13-7 

94. 5 

004 

0-298(2' 

0-553(31 

9  93(4) 

u  495(5) 

13-8 

97.5 

0035 

0'44(9i 

0-1623(41 

14  5(29) 

0-08)3) 

*  5/  V,,  IS  the  ratio  of  the  total  scattering  surface  area  to  Ihe  volume  of  ihe  sample  in  the  neutron  beam.  Ahere  the  beam  diameter  is  U  9  cm 


ing  sintering  force  (RJiincs  &  DeHoff,  19S4).  The 
isolated  pores  become  more  spherical,  achieving  radii 
which  are  greater  than  the  cross-sectional  radii  of  the 
original  channels,  yet  are  still  within  the  MSANS 
range.  As  the  number  of  isolated  pores  increases,  their 
larger  dimension  increasingly  contributes  to  the  total 
scanering,  until  it  finally  dominates  the  measured 
MSANS  intensities.  Thus,  as  seen  in  Tables  2  and  3, 
Rcfi(O)  and  /?en(<^)  begin  increasing  gradually  as  the 
transition  from  intermediate  to  final-stage  sintering 
begins.  Late  in  the  final  sintering  stage,  a  more  rapid 
increase  in  both  /?e(i(0)  and  /?e(!(°o)  is  observed  as 
the  microstructure  becomes  dominated  by  isolated 
pores. 


7.  Concluding  remarks 

The  multiple  small-angle  neui.-on-scattcring  tech¬ 
nique  was  successfully  used  to  observe  the  transition 
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Fi|.  7.  Effective  pore  radii  Tor  the  alumina  $amplei  as  a  function 
of  percent  theoretical  density  determined  from  the  MSANS 
fomalism  from  a  Porod  analysis  and  from  mercury  porosimetry 
measurements  on  the  samples  below  85%  TD,  where  open 
porosity  is  espccted. 


from  one  sintering  stage  to  another  in  the 
densification  of  alumina.  The  effective  radius  of  the 
pore-size  distribution,  as  measured  by  MSANS, 
Porod  scattering  and  mercury  porosimetry,  remained 
constant  throughout  the  intermediate  stage  of  sinter¬ 
ing,  with  a  rapid  increase  measured  by  MSANS  and 
Porod  scattering  as  densification  proceeds  beyond 
85%  of  theoretical  density.  These  results  are  consistent 
with  sintering  models  which  equate  the  intermediate 
stage  with  the  stable  decay  of  a  topological  network 
of  channels  and  junction  pores  in  the  system.  Before 
complete  densification  can  occur,  a  transition  into 
the  final  sintering  stage  takes  place  in  which  the 
porosity  forms  isolated  spherical  voids.  The  Porod 
and  MSANS  results  for  the  radii  of  the  final-stage 
voids  agree  qualitatively,  indicating  that  the  pore-size 
distribution  does  not  broaden  appreciably  as  its  peak 
advances  toward  higher  R  values  during  late-stage 
densification. 

The  agreement  among  the  MSANS,  Porod  and 
porosimetry  results  concerning  the  constant  value  of 
the  pore  radius  during  the  intermediate  stage  illus¬ 
trates  that  the  MSANS  formalism,  as  applied  to  the 
alumina  system,  is  adaptable  to  the  use  of  strings  of 
spheres  to  describe  cylindrical  pore  channels  and  is 
thus  not  limited  by  the  restriction  of  isolated  spherical 
scatterers.  The  quantitative  differences  between  the 
scattering  and  the  porosimetry  results  arise  because 
the  MSANS  and  Porod  radii  are  derived  from  a 
measure  of  all  of  the  voids  in  the  system,  with  MSANS 
radii  weighted  toward  the  larger  and  Porod  radii 
weighted  toward  the  smaller  voids.  However,  in  a 
system  consisting  of  a  homogeneous  I’Stribution  of 
pore  sizes,  si'ch  as  alumina,  the  porosimetry  radii  are 
dominated  by  the  size  of  the  small  necks  between 
pore  cha.inels. 

As  densification  proceeds  beyond  85%  TD,  the 
channels  become  pinched  off  and  then  the  transition 
from  an  open  pore  network  to  isolated  pores  occurs, 
resulting  in  an  increase  m  both  and  Rt„(oo). 

Although  even  the  largest  pores  are  expected  to  begin 
to  shrink  by  95%  TD  (Johnson,  1990),  the  scattering 
results  (Fig.  7)  indicate  that  RtniO)  and  may 
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still  be  increasing  at  97*5%  TD,  This  suggests  that 
discontinuous  grain  growth  may  be  taking  place. 
However,  the  large  uncertainties  in  the  scattering 
volume  fraction,  <p,  and  the  effective  radii,  RcrC^)  ^"d 
preclude  drawing  firm  conclusions. 

Although  the  final-stage  samples  were  originally 
intended  for  single-scattering  SANS  measurements, 
it  was  found  that  even  the  97*5%  TD  sample  was 
within  the  limit  of  applicability  of  the  MSANS  tech¬ 
nique.  There  was  significant  wavelength-dependent 
beam  broadening  and  even  the  smallest  z  value  (z  = 
40)  is  well  within  the  MSANS  formalism.  However, 
the  faa  that  <p  is  so  small  means  that  the  change  in 
r,  as  a  function  of  A  is  small  (see  Fig.  4).  This  nearly 
flat  slope  makes  the  refinement  procedure  difficult, 
resulting  in  large  uncertainties  in  the  fitted  Rei>(0) 
and  9  values.  To  improve  the  determination  of 
during  final-stage  sintering,  measurements  of  the 
higher-density  samples  are  planned  with  a  high-resol¬ 
ution  small-angle  X-ray  scattering  instrument  (Long, 
Jemian,  Weertman,  Black,  Burdette  &  Spal,  1991). 
This  instrument  can  measure  the  single-scattering 
intensities  from  pores  4  0  nm  to  10  ^.m  in  size  and 
thus  bridge  the  gap  between  those  sizes  which  can  be 
easily  measured  by  conventional  SANS  and  MSANS 
techniques. 


This  work  was  supported  in  pan  by  the  US  Army 
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Small-angle  neutron  scattering  measurements  were  used  to 
examine  the  effect  of  green  density  and  the  role  of  MgO  ad¬ 
ditive  on  the  evolution  of  the  porous  microstructure  of  alu¬ 
mina  during  intermediate-  and  final-stage  sintering.  It  was 
found  that  the  initial  connectivity  in  the  green  state  plays  a 
dominant  role  in  establishing  the  channel  diameters  during 
the  intermediate  stage  of  sintering,  and  contributes  also  to 
determining  the  onset  density  at  which  the  final  stage  of  sin¬ 
tering  begins.  The  role  of  MgO  as  a  sintering  aid  lies,  at 
least  in  part,  in  prolonging  the  stability  of  intermediate- 
stage  sintering  such  that  the  body  achieves  greater  density 
before  the  transition  to  final-stage  sintering  after  which  iso¬ 
lated  pores  are  formed.  [Key  words:  alumina,  microstruc¬ 
ture,  sintering,  characterization,  neutrons.) 

I.  Introduction 

Numerous  sintering  models  have  been  developed  over  the 
years  in  attempts  to  reliably  predict  the  kinetics  of  sinter¬ 
ing  and  to  estimate  sintering  rates.  From  a  theoretical  per¬ 
spective,  it  is  often  difficult  to  proceed  beyond  idealized 
systems  because  of  the  complexity  of  real  powder  compacts. 
Experimentally,  quantitative  characterization  is  often  ham¬ 
pered  by  Che  existence  of  a  wide  variety  of  conditions  under 
which  observations  of  the  evolving  physical  system  must  be 
made.  Although  it  has  been  difficult  to  find  a  consensus  re¬ 
garding  many  aspects  of  sintering  mechanisms,  there  is  gen¬ 
eral  agreement  that  the  driving  force  for  sintering  is  the 
surface  free  energy,  and  that  the  microstructure  characteris¬ 
tics  determine  the  specific  mechanisms  by  which  the  pore 
and  the  solid  phases  evolve  during  densification.  The  present 
investigation  was  undertaken  to  examine  the  microstructure 
evolution  during  sintering  with  the  goal  being  to  gain  insight 
into  the  mechanisms  operative  during  thermal  processing. 

Small-angle  neutron  scattering  is  a  well-developed  and  rec¬ 
ognized  technique  for  the  :,tatistically  significant  characteri¬ 
zation  of  microstructures  in  the  1  nm  to  O.I  (im  range.  This 
traditional  size  range,  which  is  too  low  to  include  most  of  the 
microstructure  relevant  to  ceramic  research,  has  recently 
been  extended  to  10  by  the  development  of  a  multiple 
neutron  scattering  formalism.  The  new  multiple  scattering 
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formalism  is  directly  applicable  to  the  important  questions 
concerning  ceramic  microstructure  evolution  during  sintering 
in  that  it  is  valid  for  samples  containing  a  large  volume  frac¬ 
tion  of  open  or  closed  pores  in  the  size  range  between  0.08 
and  10  /am.  Indeed,  the  new  technique  is  applicable  to  most, 
and  under  certain  conditions  ail,  porous  microstructures  from 
early-stage  through  intermediate-  and  final-stage  thermal 
processing. 

In  the  experiments  that  form  the  subject  of  this  paper, 
samples  of  undoped  alumina  and  of  alumina  with  0.25  wt% 
MgO  added  were  prepared  by  slip  casting,  or  alternatively  by 
cold  pressing.  For  each  green  density,  a  collection  or  suite  of 
nearly  identical  samples  were  conventionally  sintered.  During 
sintering,  these  samples  were  removed  one  by  one  from  the 
furnace  as  they  reached  densities  between  52%  and  99%  of 
theoretical,  yielding  examples  of  microstructure  from  the 
early  part  of  the  intermediate  stage  of  sintering  out  to  the 
latter  part  of  final-stage  sintering.  The  microstructure 
parameters  of  these  samples  were  derived  from  multiple  small- 
angle  neutron  scattering  measurements.  In  particular,  the 
microstructure  evolution  as  a  function  of  green  density  and  as 
a  function  of  MgO  additive  was  studied. 

II.  Multiple  Small-Angle  Neutron  Scattering  Formalism 

The  application  of  the  multiple  small-angle  neutron  scatter¬ 
ing  (MSANS)  technique'"  to  microstructure  evolution  as  a 
function  of  the  thermal  processing  of  ceramics  has  been  de¬ 
scribed  earlier.^'*  In  this  paper,  an  overview  of  the  assumptions 
involved  is  given  and  the  general  features  of  the  formalism 
are  described. 

Unlike  conventional  small-angle  neutron  scattering  (SANS), 
the  scattering  in  the  MSANS  diffraction  regime  is  dominated 
by  incoherent  multiple  neutron  scattering.  Whether  a  particu¬ 
lar  scattering  experiment  is  in  the  SANS  or  in  the  MSANS 
regime  is  determined  by  the  value  of  the  phase  shift  v  that  a 
neutron  undergoes  in  traversing  a  particle  or  pore  of  radius  R. 
The  phase  shift  depends  on  three  parameters— the  contrast 
Ap  of  the  scattering  particle  or  pore  relative  to  the  matrix, 
the  neutron  wavelength  A,  and  the  particle  size  2R: 

V  =  2ApA7?  (1) 

in  SANS,  v  1,  whereas  in  MSANS,  O.I  <  v  <  2.  A  meas¬ 
ure  of  tlie  amount  of  multiple  scattering  is  given  by  a  param¬ 
eter  I,  called  the  scattering  power,  which  is  the  sample  thick¬ 
ness,  z,  divided  by  the  neutron  mean  free  path  length,  t: 

2  -  z/f  -  1.5<^z(ApA)’/?  (2) 

d  is  the  volume  fraction  of  scattereis,  which  for  ceramic  bod¬ 
ies  more  than  50%  dense  is  I  -  d/d,,,,  where  d  is  the  sample 
density  and  d,„  is  the  theoretical  density  of  a  fully  densified 
body.  Ap  is  the  neutron  contrast  scattering-length  density*  of 
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spectra,  measured  at  each  wavelength  with  no  scatterer  in 
place,  were  also  fitted  to  en^le  a  correction*  for  instrumental 
tmjadeniog.  The  radii  of  the  pores  were  extracted  from  each 
curve  usiog  ao  empirical  formula^  for  the  description  of 
MSANS: 

riQR]  =  0.926{2  ln®“  5}®^  (for  f  a  5)  (3) 

where  f  is  given  by  Eq.  (2)  and  [. . .]  means  with  respect  to  and 
defines  the  scale  for  describing  r^.  These  values  of  R,  derived 
from  the  fitted  values  of  for  each  wavelength,  were  aver¬ 
aged  together  to  obtain  an  effective  radius  /?ctr(0)  for  the  scat¬ 
tering  pores  io  each  sample.  The  (0)  is  a  reminder  that  this 
derivation  of  effective  pore  radius  comes  from  an  analysis  of 
the  portion  of  the  scattering  curves  near  0  =*  0- 
Finally,  both  and  the  scattering  volume  fraction, 
were  reHned  simultaneously  for  each  sample  using  a  least- 
squares  fitting  method.  Because  of  the  approximate  nature  of 


Eq.  (3),  the  full  numerical  MSANS  analysis  describing  the 
dependence  of  the  radius  of  curvature  on  neutron  wave¬ 
length,  developed  in  Eqs.  2.12  to  2.1S  of  Berk  and  Hardman- 
Rhyne,'  was  used  in  the  reHnement  procedure. 

llie  measured  (symbols)  vs  wavelength  for  each  (A  the 
four  suites  of  samples  is  shown  in  Fig.  3,  along  with  theoreti¬ 
cal  values  (solid  lines)  determined  by  means  of  the  MSANS 
formalism.  The  good  ft  between  the  measurements  and  the 
theory  suggests  that  the  MSANS  formalism  is  well-suited  to 
the  analysis  of  the  porous  microstructure  of  these  samples. 

V.  Results  and  Discussion 

Results  from  the  measurements  on  the  slip-cast  alumina 
(series  13)  samples  are  considered  first.  The  Porod  results  for 
these  samples  are  given  in  Table  El,  and  the  evolution  of  the 
number  density  as  a  function  of  perceot  theoretical  density  is 


(C)  (D) 

Fig.  3.  Measured  radii  of  curvature  (symbols)  and  the  correspondine  MSANS  theoretical  (it  (solid  lines)  as  a  function  of  neutron  wave¬ 
length  for  (A)  the  t3-series.  (B)  the  V-series,  (Q  the  I5-ieries,  and  (D)  the  D-series. 
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Tabic  II.  Porod  Results  for  the  (13'Serics)  Slip-Cast  Alumina  Samples 


Sample 

Qm. 

(A'^l 

Surface  area 

Number  densiiy 
(xlO"  cm’’) 

(lim) 

*/S, 

13-2 

0.10 

15.20  i  0.2 

147  27 

0.091  ±  0.001 

3m  ±  0.04 

13-3 

0.09 

9.71  *  0.1 

66  i  15 

0.108  i  0.001 

3.60  ±  0.04 

13-4 

0.07 

5.83  *  0.03 

39.6  r  4.8 

0.108  ±  0.001 

3.60  -  0.02 

13-5 

0.065 

3.74  ±  0.02 

20.5  -  1.8 

0.120  r  0.001 

4.01  ♦  0.02 

13-1 

0.055 

2.15  ±  0.01 

7.2  ±  1.0 

0.154  2:  0.001 

5.12  ±  0.02 

13-6 

0.045 

1.218  ±  0.005 

2.84  ±  0.72 

0.185  r  0.001 

6.16  *  0.03 

13-7 

0.04 

0.553  *  0.003 

0.495  4:  0.005 

0.298  ±  0.002 

9.93  2:  0.04 

13-g 

0.035 

0.1723  r  0.0004 

0.08  ±  0.03 

0.44  2:  0.09 

14.5  *  2.9 

shown  in  Fig.  4.  The  number  of  pores  per  unit  volume  de¬ 
creases  by  several  orders  of  magnitude  as  densification  of  the 
alumina  proceeds.  Figure  5  shows  the  evolution  of  the 
MSANS-derived  pore  radius,  R,h(0)>  for  the  same  slip-cast 
alumina  samples  as  a  function  of  densification.  Although  the 
number  of  pores  is  rapidly  decreasing,  the  actual  size  of  the 
pores  is  constant  over  the  entire  range  of  densities  representa¬ 
tive  of  intermediate-stage  sintering,  from  about  55%  to  ap¬ 
proximately  87%  TD.  In  the  final  stage  of  sintering,  above 
87%  TD,  the  size  of  the  pores  increases  rapidly  despite  the 
rapid  decrease  in  their  numbers.  Also  shown  in  Fig.  5  are  the 
Porod-derived  pore  radii,  =  3<^/S„  where  (“)  indi¬ 

cates  the  large  Q  region  of  the  scattering  pattern  and  5.  is  the 
surface  scattering  area  per  unit  volume.  Finally,  results  are 
also  shown  from  MIP  measurements  over  the  range  of  densi¬ 
ties  where  the  pore  structure  was  open  enough  for  this  tech¬ 
nique  to  be  attempted.  The  /?,«(«>)  results  are  in  qualitative 
agreement  with  the  /7.n{0)  results  in  that  they  both  indicate  a 
constant  pore  radius  during  the  intermediate  stage  of  sinter¬ 
ing  and  a  rapidly  increasing  pore  size  during  the  final  stage  of 
sintering.  The  porosimetry  results  qualitatively  confirm  both 
of  the  neutron  scattering  analyses  by  again  indicating  a  con¬ 
stant  effective  pore  radius  for  the  intermediate  stage.  The  or¬ 
dering  of  sizes  from  these  three  techniques  is  as  expected 
since  the  MSANS  results  are  volume-weighted,*  emphasizing 
the  population  of  large  pores  in  the  distribution,  the  Porod 
SANS  results  are  surface-area-weighted,'  emphasizing  the 
smaller  pores  in  the  distribution,  and  the  porosimetry  results 
are  limited  by  the  narrow  necks  connecting  the  large  popula¬ 
tion  of  larger  pores. 

All  three  techniques  indicate  that  the  effective  pore  size 
and  also  the  pore  morphology  are  constant  through  the  inter¬ 
mediate  sintering  stage  of  alumina.  These  results  lend  sup^rt 
to  the  topological  decay  model  of  Rhines  and  DeHofr  in 
which  the  interconnected  network  of  pores  decays  in  a  stable 


Fig.  4.  Number  deosiiy  vs  %TD  for  the  13-series  undoped  alumina 
samplea. 


manner.  During  the  intermediate  stage,  the  interconnected 
pore  channels  decay  such  that  the  volume  fraction  to  surface 
area  per  unit  volume  (<A/S,)  remains  constant."  Consequently, 
the  pore  channels  become  fewer,  while  retaining  the  same 
diameter,  as  densiHcation  proceeds.  From  Table  II,  the  Porod 
results  for  the  ratio  of  the  volume  fraction  to  surface  area  per 
unit  volume  can  be  seen  to  have  nearly  constant  values 
throughout  the  intermediate  sintering  stage. 

Next,  MSANS  results  from  all  four  suites  of  samples  arc 
considered.  Figure  6  shows  the  evolution  of  the  pore  size  as  a 
function  of  densification  for  the  four  suites  of  samf^es,  and 
Table  III  lists  the  experimental  MSANS  parameters  and  the 
effective  pore  radii.  The  numbers  associated  with  the  left¬ 
most  symtels  in  the  figure  indicate  the  average  green  density 
of  each  suite  of  samples.  As  one  might  expect  from  the  results 
on  undoped  alumina  shown  in  Fig.  5,  the  effective  pore  radius 
was  constant  during  the  intermediate  stage  of  sintering  for 
each  of  the  suites  of  samples.  However,  the  actual  value  of 
Rrn(0)  depends  inversely  on  the  green  density  as  shown  in 
Table  IV.  The  greater  the  green  density,  the  smaller  the  effec¬ 
tive  radius  that  persists  during  intermediate-stage  sintering. 
Evidently,  once  the  size  of  the  pore  channels  is  established,  it 
is  constant  until  the  onset  of  the  final  stage  of  sintering. 
These  results  indicate  that  for  the  alumina  system  the  con¬ 
nectivity  that  is  established  in  the  green  body  determines  the 
pore  size  throughout  the  stable  intermediate  sintering  stage, 
and  once  this  size  is  determined,  it  is  unvarying  as  long  as 
isolated  porosity  has  not  yet  begun  to  form. 

At  the  onset  of  final-stage  sintering,  the  pore  radius  appar¬ 
ently  coarsens  abruptly.  Part  of  this  apparent  coarsening  may 
be  due  to  the  pinching-off  of  pore  channels  such  that  the 
larger  dimension  of  the  newly  isolated  pores  comes  increas¬ 
ingly  within  the  range  of  sizes  visible  to  the  MSANS  tech¬ 
nique.  Another  possibility  is  that  the  increase  in  the  effective 
pore  radius  may  be  smaller  than  it  seems  if  the  error  bars  on 


FIc.  5.  Effective  pore  radii  for  ihe  13-scries  alumina  as  a  function 
of  %TD  derived  from  MSANS,  from  l^odscatiering,andfroffi  MIP. 
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FI,.  6.  Effeciive  pore  radii,  /?c(i(0),  as  a  function  of  %TD  for  all 
four  series  of  samples.  The  number  appearing  near  the  leftmost 
symbols  is  the  average  green  body  density  for  each  series. 

the  measurements  could  be  reduced.  A  future  objective  of  the 
authors  is  an  investigation  of  the  final  stage  of  sintering  in 
alumina  with  improved  resolution.  The  present  objective 
which  remains  is  the  determination  of  the  influence  of  MgO 
additive  on  pore  coarsening  in  the  final  stage. 

Figure  6  indicates  that  for  the  D-series  of  samples  with 
0.2S  wt%  MgO,  the  intermediate  sintering  stage  remains 
stable  out  to  significantly  higher  density  (~93%TD)  than 
the  undoped  alumina  l3-$eries,  which  remains  stable  out  to 
— g7%TD.  Consequently,  the  pores  in  the  doped  sample  be¬ 
gin  to  coarsen  later  and  thus  do  not  suffer  the  same  degree  of 
exaggerated  growth  evident  in  the  thermal  processing  of  the 
undoped  alumina,  even  though  the  average  green  density  was 
lower  in  the  doped  sample.  Similarly,  the  transition  to  the 
final  sintering  stage  for  the  doped  sample  with  the  lowest 
green  body  density  (15-series)  occurs  no  later  in  the  densifica¬ 
tion  process  (-87%  TD)  than  the  transition  in  the  undoped  alu¬ 
mina  with  higher  average  green  density.  Evidently,  while  the 
green  density  determines  the  scale  of  the  pore  channel  size  in 
intermediate-stage  sintering  and,  in  part,  determines  when 


the  transition  to  final-stage  sintering  will  occur,  the  MgO  sin¬ 
tering  aid  can  effectively  alter  the  position  of  the  transition 
between  connected  and  isolated  piorosity.  The  transition  to 
isolated  porosity  is  delayed  in  the  doped  alumina  because  the 
additive  apparently  enhances  the  stability  of  the  intermediate 
sintering  stage  during  which  most  of  the  porosity  is  removed. 
Since  isolated  pores  begin  to  form  in  the  doped  sample  only 
after  the  remaining  porosity  is  significantly  reduced,  there  is 
less  opportunity  for  pore  coarsening  to  take  place.  Reduc¬ 
tions  in  pore  coarsening  have  also  been  observed  in  SANS 
measurements  of  plasma-sintered  MgO-doped  alumina." 

The  present  study  has  followed  the  microstructure  evolu¬ 
tion  of  porosity  as  a  function  of  sintering  to  gain  insight  into 
densification  processes.  Research  by  others'' °  has  focused  on 
the  evolution  of  the  mean  pore  spacing  or  on  the  behavior  of 
the  dihedral  angle  between  the  free  surfaces  of  grains  sharing 
a  grain  boundary.  Although  the  force  balance  is  most  directly 
related  to  the  dihedral  angle,  all  three  types  of  measure¬ 
ment-pore  size,  mean  pore  spacing,  and  dihedral  angle— are 
complementary.  In  particular,  it  was  found  that  the  mean 
pore  spacing  by  number  (see  Fig.  17  of  Ref.  12)  is  contant 
during  the  intermediate  sintering  stage  and  increases  during 
final-stage  sintering.  Thus  the  trends  observed  for  the  mean 
pore  spacing,  which  is  also  some  measure  of  the  grain  behavior, 
are  the  same  as  those  observed  in  the  pore  sizes  measured  in 
the  present  SANS  study.  Measurements*^  of  the  dihedral 
angle  in  nearly  fully  dense  alumina  and  alumina  doped 
with  0.25  wt%  MgO  found  that  while  the  median  dihedral 
angle  was  the  same  in  both  cases,  the  doped  sample  had  a 
significantly  narrower  distribution  of  dihedral  angles.  High- 
resolution  SAXS  measurements'*  currently  under  way  are 
being  used  to  examine  the  evolution  of  the  pore  size  distribu¬ 
tion,  as  a  function  of  doping,  as  alumina  approaches  full  the¬ 
oretical  density. 

VI.  Summary  and  Conclusions 

A  combination  of  small-angle  neutron  scattering  (Porod) 
measurements,  multiple  small-angle  neutron  scattering  meas¬ 
urements,  and  mercury  intrusion  porosimetry  measurements 
was  used  to  investigate  the  effect  of  green  body  density  and 
the  role  of  MgO  sintering  aid  on  the  densification  of  alumina. 
It  was  found  that  the  initial  connectivity  that  is  established  in 
the  green  body  plays  a  dominant  role  in  establishing  the  chan¬ 
nel  diameters  during  the  intermediate  stage  of  sintering.  It 


Table  III.  MSANS  Results 


Wavelength  <nm)  i  Phase  shift 


No. 

«!.TD 

Min 

Max 

/t,n  (;im) 

Min 

Max 

Min 

Max 

13-2 

54 

0.60 

0.99 

0.20  S  0.01 

152 

414 

0.13 

0.22 

13-3 

65 

0.60 

0.99 

0.19  i  0.01 

98 

289 

0.12 

0.21 

13-4 

79 

0.60 

1.09 

0.17  i  0.005 

49 

193 

0.10 

0.22 

13-5 

85 

0.60 

1.30 

0.18  £  0.01 

32 

204 

0.09 

0.28 

13-1 

89 

0.70 

1.42 

0.23  i  0.02 

47 

234 

0.17 

0.41 

13-6 

92J 

0.99 

1.42 

0J2  ♦  0.03 

96 

298 

0.35 

0.66 

13-7 

94.5 

0.99 

1.65 

0.53  i  0.05 

73 

248 

0.45 

0.93 

13-8 

97.5 

1.09 

1.65 

0.75  i  0.15 

40 

126 

0.53 

1.10 

V-1 

56 

0.64 

0.95 

0.24  ♦  0.01 

273 

519 

0.18 

0.25 

V-2 

90 

0.95 

1.26 

0J2  a:  0.01 

132 

226 

0.35 

0.45 

V-3 

% 

1.05 

U8 

0.41  -  0.02 

76 

162 

0.44 

0.72 

15-2 

60 

0.64 

0.95 

0J3  i  0.01 

261 

513 

0.26 

0.34 

15-4 

81 

0.84 

1.15 

032  i  0.01 

173 

332 

0.30 

0.42 

15-5 

87 

0.95 

1.26 

0  J7  a:  0.01 

185 

322 

0.40 

0.53 

15-6 

92 

1.05 

1.38 

0.43  i  0.02 

157 

304 

0.51 

0.75 

15-7 

95.5 

1.15 

IJl 

0.72  i  0.05 

100 

211 

0.69 

1.11 

D-1 

54 

0.53 

0.85 

0.27  ±  0.01 

206 

450 

0.17 

0.24 

D-3 

72 

0.64 

0.95 

0.25  a:  0.01 

169 

306 

0.21 

0.26 

CM 

88 

0.74 

1.05 

0.25  a:  O.Ol 

78 

161 

0.21 

0.31 

D-5 

93 

0.95 

1.26 

028  a:  0.01 

84 

152 

0.29 

0.40 

D-6 

96 

1.05 

1.51 

0J5  £  0.02 

63 

140 

0.45 

0.77 
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Tabic  IV,  Refationship  between  Green  Body  Density  and 
Intermediate-Stage  Pore  Size 


Sample  senes 

HmW  ipm) 

Green  body  density  (%) 

13- 

0.19 

53 

V- 

0.24 

51 

D- 

0.26 

50.5 

15- 

0.33 

46 

was  demonstrated  that  this  stage  is  consistent  with  the  topo¬ 
logical  model  of  Rhines  et  a/.’*'*’  It  was  also  shown  that  the 
average  green  density  plays  a  significant  role  in  determining 
the  onset  density  at  which  final-stage  sintering  begins.  The 
role  of  MgO  as  a  sintering  aid  lies,  at  least  in  part,  m  its 
ability  to  prolong  the  stability  of  the  intermediate  stage  of 
sintering  such  that  the  ceramic  body  achieves  greater  density 
before  isolated  porosity  is  formed. 

The  onset  of  final-stage  sintering  was  indicated,  in  both  the 
Porod  scattering  results  and  in  the  MSANS  results,  by  the 
appearance  of  an  increasing  size  cf  the  pore  radii  as  densifl- 
cation  proceeded.  This  upturn  in  pore  size,  which  was  not 
observed  to  turn  down  again  even  when  the  density  reached 
more  than  95%  of  theoretical,  occurred  in  the  presence  of  a 
very  significant  decrease  in  the  pore  number  density.  Since 
the  error  bars  on  the  pore  sizes  determined  by  the  present 
neutron  scattering  studies  are  too  large  to  enable  an  observa¬ 
tion  of  the  density  at  which  the  effective  pore  size  may  begin 
to  decrease,  further  investigations  of  the  final  stage  of  sinter¬ 
ing,  with  improved  experimental  resolution,  are  intended. 
Those  studies  will  attempt  to  determine  when/if  the  down¬ 
turn  in  pore  radius  occurs  as  a  function  of  densification,  and 
what  the  influence  of  MgO  additive  may  be  on  that  process. 
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Small-angle  neutron  scattering  measurements  were  used  to  examine  the  pore 
microstructure  evolution  of  glassy  silica  and  polycrystalline  alpha-alumina  as  a  function 
of  sintering.  It  was  shown  that  the  two  major  sintering  mechanisms,  viscous  flow  and 
surface  and  volume  diffusion,  lead  to  very  different  microstructure  evolution  signatures 
in  terms  of  the  average  pore  size  as  a  function  of  density.  However,  with  respect  to 
topology,  the  evolution  of  the  porosity  per  unit  surface  area  as  a  function  of  density 
is  remarkably  similar  in  the  two  systems. 


1.  INTRODUCTION 

Knowledge  of  the  microstructure  evolution  as  a 
function  of  thermal  processing  is  important  for  the  de¬ 
velopment  of  process  models  in  ceramics.  Although  it  is 
difficult  to  achieve  a  consensus  regarding  many  aspects 
of  sintering  mechanisms,  there  is  general  agreement 
that  the  driving  force  for  sintering  is  the  reduction  in 
surface  free  energy.  The  shape  of  the  minimum  free 
energy  surface  is  expected  to  influence  strongly  the 
microstructure  evolution  during  thermal  processing.  In 
effect,  the  topological  characteristics  of  a  compact  from 
moment  to  moment  influence  the  sintering  mechanisms 
by  which  the  pore  and  solid  phases  of  the  structure 
evolve.  Such  topological  characteristics  as  the  pore  sizes, 
size  distributions,  surface  areas,  and  number  distributions 
can  be  monitored  by  small-angle  neutron  scattering, 
yielding  a  direct  measure  of  the  microstnicture  that 
drives  the  sintering. 

The  present  investigation  was  undertaken  to  examine 
processing/microstructure  relationships  during  sintering 
of  two  relatively  well-understood  model  systems,  poly¬ 
crystalline  alumina,  in  which  solid-state  di^sion  mecha¬ 
nisms  dominate,  and  porous  glassy  silica,  in  which 
viscous  flow  is  dominant.  The  goal  is  to  gain  a  quan¬ 
titative  measure  of  the  structural  evolution  of  the  porous 
phase  that  takes  place  when  different  sintering  mecha¬ 
nisms  dominate  densiflcation. 

In  the  experiments  that  form  the  subject  of  this 
paper,  undoped  alumina  was  prepared  by  slip  casting  and 
microporous  silica  was  prepared  by  a  sol-gel  process. 
For  each  system,  a  collection  or  suite  of  nearly  identical 
samples  were  prepared  and  sintered.  During  sintering. 


samples  were  removed  from  the  furnace  one  by  one 
as  they  achieved  intermediate  densities  representative  of 
miaostructures  from  the  early  part  of  the  intermediate 
stage  of  sintering  to  the  late  part  of  the  final  stage  of 
sintering.  Microstructure  parameters  of  these  samples 
were  derived  from  small-angle  neutron  scattering  and 
multiple-small-angle  neutron  scattering,  as  required  by 
the  sizes  and  volume  fraaions  of  die  pores.  This  process 
yielded  a  quantitative  signature  of  the  structural  evolu¬ 
tion  of  crystalline  and  glassy  materials  as  a  funaion  of 
sintering. 

II.  SMALL-ANGLE  NEUTRON  SCATTERING 

Small-angle  neutron  scattering  (SANS)  is  a  well- 
recognized  technique  for  the  statistically  significant  char¬ 
acterization  of  microstructure  in  the  1  nm  to  0.1  fim  size 
range.'  This  traditional  size  range,  which  is  too  low  to 
include  many  of  the  microstructural  features  relevant  to 
ceramic  research,  has  recently  been  extended  to  10  pm 
by  the  development  of  a  multiple-small-angle  neutron 
scattering  (MSANS)  formalism.^  The  new  multiple- 
scattering  formalism  is  valid  for  samples  containing  a 
large  volume  fraction  of  open  or  closed  pores  in  the  size 
range  between  0.08  and  10  /rm.  For  the  microporous 
silica  used  in  this  study,  MSANS  was  applicable  to 
samples  in  the  intermediate  stage  of  sintering.  For  the 
alumina  used  in  this  study,  MSANS  measurements  were 
valid  from  early  intermediate  stage  through  to  the  final 
stage  of  thermal  processing.  The  reason  that  MSANS 
could  not  be  used  for  all  of  the  structures  in  this 
study  is  that  an  important  change  takes  place  in  the 
microstnicture  morphology  within  one  of  the  systems. 
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The  application  of  the  MSANS  technique  to  micro- 
structure  evolution  in  ceramics  as  a  function  of  sintering 
has  been  described  earlier/-^  In  this  paper,  the  earlier 
measurements  are  reviewed  and  extended,  and  conclu¬ 
sions  are  drawn  based  on  the  existence  of  either  glass- 
and-void  or  crystalline-and-void  phases  in  the  body. 

Unlike  conventional  SANS,  the  scattering  in  the 
MSANS  diffraction  regime  is  dominated  by  incoherent- 
multiple-neutron  scattering.  Since  the  radii  of  the  scat¬ 
tering  pores  are  between  0.08  and  10  ^m,  which  is  of  the 
same  order  of  magnitude  as  the  neutron  mean  free  path 
length,  the  neutrons  scatter  ffom  only  one  pore  at  a  time, 
even  though  they  may  scatter  many  times  before  leaving 
the  sample.  This  means  that  there  is  no  contribution  to 
the  measured  scattering  from  interference  between  pores. 

Whether  a  particular  scattering  experiment  is  in  the 
SANS  or  in  the  MSANS  regime  is  determined  by  the 
phase  shift  u  that  a  neutron  undergoes  in  traversing  a  par¬ 
ticle  of  radius  R.  The  phase  shift  depends  on  three  param¬ 
eters:  the  contrast^  Ap  of  the  pore  relative  to  the  matrix, 
the  neutron  wavelength  A,  and  the  particle  size  2R. 

i/  =  2ApAi2.  (1) 

In  SANS,  <  1,  whereas  in  MSANS,  0.1  ^  i/  <  2. 

A  measure  of  the  amount  of  multiple  scattering  is 
given  by  the  parameter  z,  called  the  scattering  power,^ 
which  is  the  sample  thickness  z  divided  by  the  neutron 
mean  free  path  length  I 

z  =  zfl  —  \.o4>z{^pX)^ R .  (2) 

4>  is  the  volume  fraction  of  scatterers,  which  for  ceramic 
bodies  more  than  50%  dense  is  (1  -  d/d^\  where  d 
is  the  sample  density  and  d,h  is  the  theoretical  density 
of  a  fully  densified  body.  In  the  present  case,  d>  is  the 
porosity.  Ap  is  the  neutron  scattering  length  density^  of 
the  pores  relative  to  the  matrix  (that  is,  the  particles) 
mentioned  above,  A  is  the  neutron  wavelength,  and  R  is 
the  mean  effective  radius  of  the  pores.  From  a  materials 
perspective,  the  greater  the  number  of  pores  and  the 
larger  their  size,  the  greater  the  multiple  scattering. 
However,  even  if  the  pore  volume  fraction  is  small,  as  in 
the  final  stage  of  sintering,  if  the  remaining  pores  have 
large  radii,  multiple  scattering  measurements  can  still 
be  carried  out  successfully,  provided  that  the  samples 
are  thick  enough  and  long  wavelength  neutrons  (of  the 
order  of  1.2  to  2  nm)  are  available.  For  sintering  of 
microporous  silica,  the  pore  structure  within  the  samples 
representative  of  the  final  stage  could  not  be  measured 
by  MSANS,  and  the  structural  evolution  measurements 
had  to  be  completed  using  SANS  measurements.  The 
absence  of  MSANS  for  the  final  stage  densification  of 
silica  implies  that  there  are  no  large  spherical  pores  left, 
as  will  be  explained  below.  For  conventional  sintering  of 
alumina,  the  entire  pore  evolution  from  early  intermedi¬ 
ate  stage  to  late  final  stage  can  be  observed  by  means  of 


MSANS  because  even  the  residual  pores  near  the  end  of 
final-stage  sintering  can  be  approximated  by  spheres  and 
are  large  enough  to  be  measured.  It  was  found  that  the 
measurement  techniques  that  can  be  used  yield  important 
clues  as  to  the  differences  that  will  be  observed  between 
the  sintering  of  these  crystalline  and  glassy  systems. 

The  MSANS  analysis  assumes  that  the  scatterers  can 
be  described  as  spheres  with  no  internal  structure  in  a 
featureless  matrix.  Of  course,  the  pores  in  both  systems 
have  no  internal  structure.  For  the  silica  case,  the  single 
colloidal  particles  and  the  particle  agglomerates  have 
been  shown*  to  be  spherical.  The  voids  between  the  par¬ 
ticles  for  porosities  between  85%  and  50%  were  shown 
in  Ref.  6  to  be  reasonably  well  approximated  by  spheres. 
The  point  at  which  the  pores  may  become  isolated  could 
not  be  determined  in  that  work.  However,  the  matrix  is 
always  glassy,  and  thus  it  is  relatively  featureless.  For 
alumina,  the  crystalline  matrix  contains  features  such  as 
dislocations  which  will  cause  scattering.  Such  structures 
do  not  generally  fall  into  the  same  size  range  as  the  pores 
and  are  thus  outside  the  region  of  sensitivity  of  MSANS. 
For  alumina,  the  spherical  approximation  for  the  shape 
of  the  isolated  pores  is  applicable  to  the  final  stage  of 
sintering,  but  there  is  some  question  as  to  how  well  the 
method  will  perform  during  the  intermediate  stage  of 
sintering  when  the  pores  form  an  interconnected  network 
of  cylindrical  channels.  Relevant  to  this  concern,  it  was 
found'*  that  the  spherical  approximation  is  valid  provided 
that  the  cylindrical  channels  can  be  modeled  by  strings 
of  spheres.  This  will  be  true  when  the  length  dimension 
of  the  cylinders  is  beyond  the  range  of  sizes  measured 
by  MSANS.  For  intermediate  stage  sintering  of  alumina, 
the  length  of  the  interconnected  pore  channels  may 
be  well  beyond  10  pm.  Supporting  evidence  for  the 
validity  of  the  strings-of-spheres  approximation  for  alu¬ 
mina  comes  from  mercury  intrusion  porosimetiy  (MIP) 
measurements  on  the  same  intermediate  stage  samples. 
In  MIP,  a  cylindrical  morphology  is  assumed  for  the 
pore  shape.  Since  the  MIP  results*  and  the  MSANS 
results  showed  the  same  behavior,  this  can  be  taken  as 
evidence  that  each  approximation  serves  appropriately 
in  the  respective  method. 

In  conventional  SANS,  the  width  of  the  scatter¬ 
ing  curve  as  a  function  of  scattering  wave  vector  Q 
(where  |Q|  =  (47r/A)  sin  and  is  the  scattering 
angle)  is  independent  of  the  incident  neutron  wavelength. 
In  MSANS,  the  width  of  the  scattering  curve  changes  as 
a  function  of  wavelength;  that  is,  there  is  wavelength- 
dependent  beam-broadenmg.  Thus,  with  the  phase  shift 
between  0.1  and  2,  with  the  pore  size  between  0.08  and 
10  pm,  and  with  z  ^  10,  the  wavelength  dependence  of 
the  MSANS  data  can  be  used  to  obtain  use^l  pore  size 
information.  By  studying  the  samples  over  a  range  of 
incident  neutron  wavelengths,  the  statistics  of  the  data 
refinement  are  significantly  improved.  The  parameters 
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required  for  the  data  analysis  are  sample  thickness  z, 
material  contrast  Ap,  the  scattering  volume  fraction 
(which  is  the  same  as  the  porosity),  and  the  wavelength 
dependence  of  the  beam  broadening  over  at  least  three, 
and  preferably  more,  wavelengths. 

Even  under  scattering  conditions  where  the  low- 
wave-vector  portion  of  the  scattering  curve  near  0  =  0  is 
dominated  by  multiple  scattering,  the  large-wave-vector 
portion  of  the  scattering  curve,  defined  by  QR  >  10, 
follows^  single  particle  Porod^  scattering.  This  means 
that  the  usual  analysis  of  the  large-^  portion  of  the 
scattering  curve  can  be  carried  out  independent  of  the 
existence  of  MSANS  at  low-^.  This  is  very  useful, 
because  a  Porod  analysis  yields  values  for  the  total 
surface  scattering  area  of  the  pores. 

dft  VQ^  '  ^ 

V  is  the  sample  volume  in  the  neutron  beam  and  5  is  the 
total  surface  area  of  the  scattering  pores.  From  this  the 
Porod  constant  P  can  be  defined  as 

P  =  2i:{ApfiS/V),  (4) 

and  one  has  a  direct  measure  of  the  total  surface  area  of 
the  pores  per  unit  volume  that  is  independent  of  model 
for  the  pore  shape.  If  the  pores  can  be  assumed  to  be 
spherical,  one  has  another  measure  of  their  size  since 
R  =  30(V7S),  and  one  can  also  obtain  a  value  for  their 
number  density  from  N/>  =  (5/V)(l/(4jr/?^)]. 

III.  SAMPLES 

The  silica  samples  were  prepared  by  the  Shoup* 
technique  from  85  wt.  %  potassium  silicate  and  15  wt.% 
colloidal  silica.  Detailed  information  on  sample  prepara¬ 
tion  is  given  elsewhere.’  The  porosity  of  the  nominally 
alkali-free  (<I(X)  parts  in  10*  Na  and  K)  precursor  body 
was  84%.  ^ch  of  the  samples  was  subjected  to  a  firing 
temperature  above  the  onset  of  sintering,  which  occurs 
between  1420  and  1470  K,  depending  on  the  average 
pore  size  of  the  dried  precursors.  Crystallization  was 
carefully  avoided  in  the  sintering  of  this  material.  The 
densities  of  the  small-angle  neutron  scattering  samples 
were  between  60%  and  98%  of  theoretical  density  (TD), 
where  the  TD  for  silica  is  2.20  g  cm'^.  Thus  the  samples 
are  representative  of  intermediate-stage  to  final-stage 
thermal  processing.'*^  The  silica  sample  characteristics 
are  given  in  Table  I. 

The  alumina  samples  were  prepared  from  Baikowski 
alumina  powder  (CR6,  6  m^  g'‘),  and  were  slipcast  to 
an  average  green  density  of  52.8  ±  0.4%  TD,  where 
TD  for  alumina  is  3.965  g  cm'^.  Detailed  information  on 
sample  preparation  is  given  elsewhere.^  During  sintering, 
the  samples  were  removed  from  the  furnace  one  at  a 
time,  yielding  examples  of  microstructure  between  56 


TABLE  I.  Silica  samples. 


Sample 

number 

Porosity 

<5 

Thickness. 

cm 

Densily, 
g  cm'^ 

1 

0.40 

0.487 

1.30 

2 

0J5 

0.142 

1.42 

4 

0.30 

0.123 

1.56 

5 

0.26 

0.122 

1.63 

6 

0.05 

0.486 

2.09 

7 

0.05 

0.113 

2.10 

9 

om 

0.130 

2.15 

and  98%  TD.  The  alumina  sample  characteristics  are 
given  in  Table  11. 

IV.  MULTIPLE  SCATTERING  MEASUREMENTS 

The  neutron  scattering  measurements  were  per¬ 
formed  at  the  20  MW  research  reactor  at  the  National 
Institute  of  Standards  and  Technology  (NIST).  The  NIST 
SANS  facility"  makes  use  of  a  velocity  selector  for 
choosing  the  mean  wavelength  of  the  incident  neutrons, 
which  have  a  wavelength  spread  of  AA/A  =  0.25. 
Copious  neutron  fluxes  are  available  in  the  wavelength 
region  between  0.5  and  2.0  nm  due  to  a  cryogenic 
moderator  installed  in  the  reactor  core.  The  neutron  beam 
was  collimated  with  a  12  mm  aperture  following  the 
velocity  selector  and  an  8  mm  aperture  preceding  the 
sample  position.  Scattered  neutrons  were  detected  by 
a  64  X  64  cm*  position-sensitive  detector  divided  into 
128  X  128  pixels.  The  data  were  circularly  averaged  to 
produce  one-dimensional  intensity  liQ)  versus  scattering 
vector  Q  curves. 

In  conventional  SANS,  a  beamstop  is  used  in  the 
center  of  the  scattering  pattern  to  protect  the  detector 
from  the  unscattered  transmitted  beam.  This  transmitted 
beam  is  10^  to  10’  times  greater  in  intensity  than  the  scat¬ 
tered  beam.  For  MSANS,  measurements  were  performed 
with  no  beamstop  in  place  since  the  intense  transmitted 
beam  associated  with  SANS  is  absent.  Performing  the 
measurements  without  a  beamstop  ensures  that  the  beam 


TABLE  tl.  Alumina  samples. 


Sample 

number 

Potosiiy 

<P 

Thickness, 

cm 

Density, 
g  cm'^ 

t3-2 

0.46 

1.019 

2.15 

13-3 

0J5 

0.988 

2.59 

13-4 

0.21 

0.894 

3.t5 

t3-5 

0.t5 

0.9t4 

3J9 

13-1 

0.11 

0.851 

3.55 

t3-6 

0.075 

0.852 

3.69 

t3-7 

0.055 

0  840 

3.77 

t8-8 

0.025 

0.804 

3.89 

2T06 
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center  can  be  located  with  good  precision,  and  enables 
an  accurate  determination  of  the  beam-broadening  of  the 
I{Q)  curve  at  Q  =s  0.  Measurement  times  were  of  the 
order  of  10  min  to  1  h,  depending  on  the  amount  of 
multiple  scattering  present. 

MSANS  measurements  were  made  at  a  minimum  of 
five  wavelengths  per  alumina  sample.  Multiple  scattering 
curves  from  one  of  the  alumina  samples  (13-5)  are  shown 
in  Fig.  1,  where  all  of  the  curves  have  been  normalized  to 
/(O)  =  1.  The  feature  at  0  ~  0.05  nm'*,  which  appears 
in  all  of  the  broadened  spectra,  is  due  to  a  change 
in  the  detector  sensitivity  near  2  =  0.  To  avoid  this 
region,  Gaussian  fits  to  the  data  were  made  beginning  at 
85-90%  /(O)  and  extending  out  to  approximately  45% 
1(f)).  Below  40%  /(O),  as  the  theory  predicts,  the  curves 
are  not  well-described  by  a  Gaussian.  For  the  silica 
samples,  the  region  that  was  distorted  by  the  changes  in 
detector  sensitivity  was  larger  because  a  different  portion 
of  the  detector  was  used.  Therefore,  the  silica  curves 
could  be  fit  only  between  75%  and  40%  of  1(0). 

The  wavelength  dependent  beam-broadening  of  the 
MSANS  curves  is  described  as  a  radius  of  curvature, 
Tc,  at  2  =  0  jn  the  scattering  curve,  in  keeping  with 
the  earlier  literature.^  For  each  scattering  curve,  a 
radius  of  curvature  was  extracted  from  the  fitted  standard 
deviation  a,  such  that  =  (l/la)^^^.  Blank  spectra, 
measured  at  each  wavelength  with  no  sample  in  place, 
were  also  measured  to  enable  a  conection^^  for  instru¬ 
mental  broadening.  The  radii  of  the  pores  were  extracted 
from  each  curve  using  an  empirical  formula^  for  the 
description  of  the  MSANS: 

r,.[2H]  =  0.926{z  ln°  2}“  ®  for  2  >  10 ,  (5) 


where  2  is  given  by  Eq.  (2)  and  (...)  means  “with 
respect  to”  and  defines  the  scale  for  describing  r^.  The 
values  of  R,  derived  from  the  fitted  values  of  r^,  are 
averaged  together  to  obtain  an  effective  radius  /?efr(0) 
for  the  scattering  pores  in  the  sample.  The  (0)  is  a 
reminder  that  this  value  of  the  pore  radius  comes  from 
the  2  =  0  portion  of  the  scattering  curve.  Later  in  this 
paper,  another  pore  radius  will  be  derived  from  the  large 
2  (Porod)  portion  of  the  scattering  curves,  and  it  will  be 
referred  to  as  Rettfoc). 

To  complete  the  MSANS  analysis,  both  ^eir(0)  and 
the  scattering  volume  fraction  tP  are  refined  simulta¬ 
neously  for  each  sample  using  a  least-squares  fitting 
procedure.  Due  to  the  approximate  nature  of  Eq.  (5),  the 
full  numerical  MSANS  analysis  describing  the  depen¬ 
dence  of  the  radius  of  curvature  on  neutron  wavelength, 
developed  in  Eqs,  (2.12)-(2.15)  of  Ref.  2,  was  used  in 
the  refinement  procedure. 

The  measured  r^’s  as  a  function  of  wavelength  for 
the  silica  samples  are  shown  in  Fig.  2,  and  the  measured 
re’s  versus  A  for  the  alumina  samples  are  shown  in 
Fig.  3.  The  lines  in  Figs.  2  and  3  arc  the  theoretical 
curves  determined  by  means  of  the  MSANS  formalism. 
The  data  for  silica  sample  4,  with  <(>  =  0.3  was  too  close 
to  the  data  for  silica  sample  5  with  4>  =  0.26,  so  it 
was  omitted  from  the  figure  for  clarity.  Nevertheless, 
both  figures  indicate  good  agreement  between  the  theory 
and  the  measurements,  demonstrating  that  the  MSANS 
formalism  is  well-suited  to  the  analysis  of  these  porous 
microstructures. 

It  is  interesting  to  note  that,  regardless  of  whether 
the  phase  of  the  material  is  glassy  or  crystalline,  the 


FIG.  2.  Radii  of  curvature  for  the  silica  samples  derived  from  the  data 
(symbols)  and  calculated  from  the  theory  (doited  lines)  as  a  function 
of  neutron  wavelength. 
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FIG.  3.  Radii  of  curvature  for  the  alumina  samples  derived  from 
the  data  (symbols)  and  calculated  from  the  theory  (solid  lines)  as 
a  function  of  neutron  wavelength. 


versus  A  curves  have  the  steepest  slopes  for  the  data  from 
the  samples  with  the  largest  volume  fraction  of  pores, 
and  the  shallowest  slopes  for  the  data  from  the  samples 
with  the  least  porosity.  This  trend  was  observed  also'^ 
in  MSANS  data  from  alumina  samples  that  had  been 
doped  with  0.25  wt.  %  MgO.  The  results  indicate  that 
while  the  size  of  the  pores  may  be  increasing,  decreasing, 
or  remaining  constant  as  the  porosity  is  decreasing,  the 
general  trend  of  the  versus  A  curves  is  observed  to 
behave  as  shown  in  Figs.  2  and  3.  Thus,  even  though 
both  and  /7,cK0)  arc  implicit  in  these  curves,  the  slope 
is  dominated  by  the  volume  fraaion  of  pores  <t>. 


V.  RESULTS  AND  DISCUSSION 

Results  from  the  sintering  of  microporous  silica 
are  considered  first.  The  Porod  results  for  silica  are 
given  in  Table  111.  The  Porod  analysis  yields  directly 
a  value  for  the  total  surface  scattering  area.  This 
surface  area,  normalized  to  unit  volume  of  sample  in 
the  neutron  beam,  is  generally  seen  to  decrease  as 
densification  proceeds.  Assuming  a  spherical  shape  for 
the  pores,  a  value  for  their  average  radius  f?eff(oc) 
can  be  derived,  and  this  is  also  given  in  Table  111. 
The  average  radius  is  first  seen  to  coarsen  and  then 
to  shrink.  Finally  the  radius  appears  to  become  very 
large,  which  is  hard  to  understand  unless  this  is  taken 
as  evidence  that  the  spherical  approximation  fails  for 
the  pore  shape  in  silica  over  95%  TD.  Since  it  will 
be  demonstrated  that  the  pore  morphology  in  final- 
stage  silica  is  disk-like,  the  last  two  numbers  in 
Table  111  are  not  physically  meaningful  and  are  left  for 
reference  purposes  only.  Another  parameter,  the  pore 
volume  fraction  divided  by  the  surface  area  normalized 
to  unit  volume,  is  also  shown  in  Table  III.  Since 
this  is  a  topological  parameter,  all  of  the  values  for 
0(S/V)  shown  in  the  table  are  physically  relevant. 
More  will  be  said  about  the  meaning  of  this  param¬ 
eter  below. 

Near  Q  =  0,  MSANS  was  used  to  measure  the 
average  pore  sizes  of  microporous  silica.  During  the 
intermediate  stages  of  thermal  processing  there  was  a 
considerable  amount  of  wavelength-dependent  broaden¬ 
ing;  however,  there  was  no  detectable  broadening  in  the 
late  stage  samples,  even  with  the  longest  wavelength 
neutrons  (2  nm)  and  the  thickest  samples.  In  an  earlier 
study.'*  it  was  demonstrated  that  pores  with  radii  between 
0.08  and  3  pm  would  have  produced  detectable  broad¬ 
ening  using  the  longest  wavelength  neutrons  available. 
The  absence  of  such  broadening  in  the  final  stage  indi¬ 
cates  that  the  large  sphere-like  voids  present  during  the 
intermediate  stage  of  sintering  were  completely  removed 
before  the  silica  reached  95%  TD.  To  observe  the  pore 
population  in  the  late  stage  of  sintering,  measurements 
were  made  by  means  of  SANS.  In  the  SANS  analysis, 
the  scattering  pore  size  was  estimated  by  fitting  the  small 


TABLE  in.  Silica  Porod  results. 


Sample 

number 

Porod  consianu 

X  10“  cm*’ 

S/V 

X  10*  cm*’ 

Porosity 
d>.  % 

H(S/V) 

X  1(3^  cm 

Din 

1 

117  ±  0.34 

16.98  ±0.40 

40 

2.36  ±  0.06 

71  ±  2 

2 

8.42  ±0.14 

11.26  ±0.18 

35 

3.11  ±  0.05 

93  ±  2 

4 

9J4  ±0.12 

12.76  ±0.16 

30 

2  J5  ±  0.03 

70  ±  1 

5 

7.88  ±0.04 

1034  ±0.05 

26 

2.46  ±  0.01 

74  ±  03 

6 

0.283  ±  0.002 

0.378  ±  0.003 

5 

13.2  ±  1.0 

397  ±  28 

9 

0.118  ±  0.001 

0.156  ±0.001 

2 

12.8  ±0.8 

384  ±32 
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Q  scattering  data  to  the  Gutnier  equation'"* 

(6) 

where  Rg  is  the  radius  of  gyration  of  the  pores  and  /(O)  is 
the  extrapolated  scattering  intensity  at  0  =  0.  The  SANS 
results  on  the  late  stage  samples  revealed  the  presence 
of  small  pore  populations  with  radii  of  gyration  33  nm 
and  less.  The  MSANS  and  SANS  results  are  shown  in 
Table  IV.  It  can  be  seen  that  the  pores  first  coarsen 
and  then  shrink  abruptly  as  the  silica  densifies.  These 
results  support  a  model’  proposed  for  this  system  in 
which  there  are  two  distinct  populations  of  pores.  There 
is  a  population  of  small  pores  that  exist  within  loosely 
aggregated  clusters  or  “domains”  of  closely-packed  fine 
silica  particles.  These  pores  within  domains  are  actually 
voids  between  touching  particles.  Also  there  is  a  popula¬ 
tion  of  large  pores  that  are  found  between  domains.  The 
latter  population  of  pores  dominates  the  MSANS  data. 
The  small  pores  are  less  stable  and  thus  are  eliminated 
first,  causing  the  domains  to  densify  and  the  interdomain 
pores  to  coarsen.  When  the  intradomain  pores  have 
disappeared,  the  interdomain  pores  can  stop  coarsening 
and  begin  to  shrink.  These  small-angle  neutron  scattering 
results  indicate  that  it  is  possible  to  determine  exactly 
when  domain  densification  is  complete. 

Comparing  the  Porod  and  the  MSANS  and  SANS 
results  for  pore  radius,  one  can  see  that  the  Porod- 
derived  radius,  which  is  a  measure  of  all  of  the  porosity 
but  weighted  toward  the  smaller  pore  population,  begins 
to  decrease  earlier  than  the  MSANS-derived  radius. 
This  can  be  taken  as  evidence  of  the  existence  of  the 
population  of  small  pores  that  is  beyond  the  range  of 
visibility  of  the  MSANS.  With  regard  to  the  final-stage 
data,  the  lack  of  MSANS  indicates  that  there  are  no  large 
pores  left  while  the  Porod  results  imply  the  presence 
of  large  pores.  The  apparent  contradiction  is  removed 
when  one  considers  the  question  of  the  final  stage  silica 
pore  shape.  For  the  final  stage,  consistency  among  all  the 
results  is  achieved  only  when  the  pore  shape  is  taken  as 
disk-like.  The  scattered  intensity  from  disk-like  pores  of 


/(g)  =  7(0)  exp  I 


finite  thickness,  T.  and  a  much  larger  diameter  can  be 
written  as*'’ 

=  (7) 

where  A  is  a  constant  and  /,(Q)  is  the  “thickness  factor”. 
Analagous  to  Eq.  (6),  /,(Q)  can  be  approximated  at  small 
Q  by 

/,(g)  =  4(0)exp(-g2fln  (8) 

where  R,  =  T/  \^T2.  Thus,  by  fitting  I{Q)  •  Q'  vs  Q- 
at  small  Q,  T  can  be  obtained.  The  results  are  given 
in  Table  IV.  Thus,  if  the  final-stage  porosity  in  the 
microporous  silica  is  taken  to  be  disk-like,  the  radii 
of  gyration  measured  by  means  of  SANS  indicate  disk 
thici^esses  of  S50  nm. 

Finally,  electron  microscope  micrographs  were  taken 
from  fracture  surfaces  of  the  74%  and  the  95%  dense 
silica.  These  are  shown  in  Figs.  4(a)  and  4(b).  At  74%, 
the  porosity  is  still  interconneaed  by  large  cylindrical 
pores,  but  some  narrow  porosity  has  also  formed.  By 
95%  dense,  the  microstructure  contains  no  large  pores 
and  only  the  narrow  collapsed  porosity  is  seen.  These 
micrographs  confirm  the  change  in  pore  morphology 
observed  by  MSANS  and  SANS. 

Next  the  alumina  results  are  considered.  The  Porod 
results  for  alumina  are  given  in  Table  V.  The  surface  area 
per  unit  volume  decreases  monotonically  as  densification 
proceeds.  However,  the  pore  radius,  assuming  spheres, 
is  seen  to  be  relatively  constant  during  the  intermediate 
stage  of  sintering,  and  increases  rapidly  during  the  final 
sintering  stage. 

The  MSANS  results  for  alumina  are  given  in 
Table  VI,  where,  unlike  silica,  MSANS  characterization 
could  be  used  throughout.  It  can  be  seen  from  the 
table  that  the  effective  size  of  the  pores  is  almost 
constant  over  the  entire  range  of  densities  representative 
of  intermediate  stage  sintering,  from  approximately  55% 
TD  to  about  87%  TD.  Above  87%  TD,  the  sizes  of 
the  pores  increase  rapidly,  despite  the  sharp  decrease 
in  their  numbers.  This  trend,  i.e.,  constant  pore  radius  in 


TABLE  IV.  MSANS  and  SANS  results  from  silica. 


Sample 

number 

Porosity 

(1>.  % 

/im 

T 

nm 

j 

1/ 

min 

max 

min 

m9x 

1 

40 

0.176  ±  0.01 

. . . 

32 

72 

0.11 

0.14 

2 

35 

0.225  ±  0.015 

21 

38 

0.19 

0.34 

4 

30 

0.223  ±  0.015 

22 

49 

0.21 

0.36 

5 

26 

0.118  ±  0.02 

16 

31 

0.15 

0.28 

6 

5 

<0.08 

7 

5 

... 

45  ±  6 

9 

2 

<0.08 

49  ±  4 
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(C) 

FIG.  4.  ScauiivdsctreomicrDgnphsof  (ai  a  fraciurv  surfuce  uf  74'<  dvn>».-  microporous  silica,  lb)  a  fracture  surface  of  95%  dease  mKCoporous 
iOka,  and  (c)  a  poHahed  sutfece  of  95'X  dense  alumina. 


die  intomediate  sintering  stage  and  rapidly  increasing 
pore  redins  in  the  final  sintering  stage,  is  alsti  seen  ior 
Rta(jxi)  whidt  was  derived  front  Ponxl  scattering,  w  here 
I  Spherical  approximation  was  made  lor  the  ptirc  shape. 
The  qualitative  agreement  between  the  Porod  and  the 
MSA]^  results  for  average  pore  radius  suggests  thai. 
uiiinra  silica,  diere  is  only  one  pore-si/;e  disirihuiion 


present  in  alumina.  The  MSANS  results  are  weighted 
toward  the  larger  pores  in  the  distributioo  and  the 
Porod  results  are  weighted  toward  die  imaller  pores  in 
;ne  distribution,  but  it  appears  that  diere  is  only  one 
distribution  present. 

The  constant  pore  radius  in  the  intermediate  sin¬ 
tering  stage  of  alumina  lends  siqiport  to  tbe  topolog- 
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TABLE  V.  Alumina 

Porod  results. 

Sample 

number 

Porod  constant 

X  lO-"  cm-‘ 

5/V 

X  10"'  cm"' 

Porosity 
<p.  '7.' 

t>l(SIV) 

X  lO"*  cm 

nm 

13-: 

31.3  i  0.4 

15.2 

±  0.2 

46 

o 

p 

91  *  1 

13-3 

20.0  i  0.2 

9.71 

±  0.1 

35 

3.60  ±  0.04 

108  i  1 

13-4 

12.0  =  0.07 

5.83 

±  0  03 

21 

3.60  ±  0.02 

108  ±  1 

13-5 

7.70  ±  0.04 

3.74 

±  0.02 

15 

4.01  ±  0.02 

120  ±  1 

13-1 

4.42  ±0j03 

2.15 

±  0.01 

11 

5.12  ±  0.02 

154  ±  1 

13-6 

2.51  ±0.011 

1.218 

±  0.005 

7.5 

6.16  ±  0.03 

185  ±  1 

13-7 

1.14  ±  0  006 

0.553 

±  0.003 

5.5 

9.93  ±  0  04 

298  ±  2 

13-8 

0.355  ±  0.0009 

0.5123  ±  0.0009 

2.5 

14i  ±  0.9 

440  ±  90 

ical  decay  model  of  Rhines  and  DeHoff'*  in  which 
the  interconnected  channel  network  of  pores  decays 
in  a  stable  manner.  During  the  intermediate  stage,  the 
interconnected  pore  channels  are  expected  to  decay 
such  that  the  volume  fraction  to  surface  area  per  unit 
volume,  4>1{S/V),  remains  constant,”  as  can  be  seen 
from  Table  V.  Consequently,  the  pore  channels  become 
fewer,  while  retaining  the  same  diameter,  as  densiheation 
proceeds.  At  the  onset  of  final-stage  sintering,  the  pore 
radius  coarsens  abruptly.  Part  of  this  apparent  coarsening 
may  be  due  to  the  pinching  off  of  pore  channels  such 
that  the  larger  dimension  of  the  newly-isolated  pores 
comes  increasingly  within  the  range  of  sizes  visible  to 
the  MSANS  technique.  The  error  bars  on  these  late-stage 
radii  are  large  because  it  is  difricult  to  fit  the  shallow 
slopes  of  the  relevant  versus  A  curves.  The  authors 
have  recently  attempted  to  improve  the  precision  of 
these  results  with  independent  measurements  using  high- 
resolution  small-angle  x-ray  measurements on  the  same 
late-stage  alumina  samples.  The  new  results  confirm 
the  average  pore  sizes  and  the  trends  shown  here.  A 
micrograph  of  a  polished  surface  of  the  95%  dense 
sample  is  shown  in  Fig.  4(c). 

The  behavior  of  the  total  surface  area  of  the  pores, 
as  a  function  of  densification,  is  shown  in  Fig.  5  for  both 
silica  and  alumina.  Although  there  are  some  quantitative 
differences  below  80%  TD,  it  is  interesting  to  see  that  the 
total  surface  area  of  the  pores  in  the  glassy  system  and 


in  the  crystalline  system  falls  on  the  same  densification 
curve.  Next,  the  parameter  o/iSfV),  which  was  taken 
above  to  be  an  indicator  of  the  existence  of  a  topologi¬ 
cally  stable  sintering  stage  in  intermediate  stage  sintering 
of  alumina,  is  compared  for  the  uensification  of  alumina 
and  silica  in  Fig.  6.  The  topology  of  the  microstructure 
evolution  of  silica  and  alumina  are  on  the  same  general 
curve.  This  is  a  remarkable  result.  The  pore  morphology 
for  the  polycrystalline  alumina  goes  from  interconnected 
cylindrical  channels  to  pinched-off  sphere-like  shapes; 
the  pore  morphology  for  the  glassy  silica  goes  from 
interconnected  cylindrical  channels  to  a  collapsed  disk¬ 
like  shape  that  is  largely  interconnected  even  at  95% 
dense.  Nevertheless,  the  topological  parameter  <i/(s/v) 
evolves  in  the  same  general  manner  during  the  sintering 
of  both. 

The  SANS  and  MSANS  results  for  the  pore  radius 
as  a  function  of  densification  are  compared  for  the  two 
systems  in  Fig.  7.  In  this  figure,  the  differences  in  signa¬ 
tures  of  the  sintering  of  a  glass  undergoing  viscous  flow 
and  crystalline  material  densifying  mainly  by  surface  and 
volume  diffusion  are  evident.  The  glassy  silica  flows 
such  that  the  smaller  voids  are  first  removed,  causing 
the  temporary  coarsening  of  the  largest  pores.  The  point 
at  which  all  the  intradomain  pores  have  vanished  can 
easily  be  seen.  After  that  point,  the  large  pores  shrink 
rapidly  and  eventually  collapse  as  the  silica  sinters  to 
full  theoretical  density.  The  alumina  densifies  by  the 


TABLE  VI.  MSANS  results  for  alumina. 


Sample 

number 

Porosity 
<t>.  % 

ficiriO) 

e 

u 

min 

mu 

min 

max 

13-2 

46 

0.20  ±  0.01 

152 

415 

0.13 

0.22 

13-3 

35 

0.19  ±  0.01 

98 

289 

0  12 

0.21 

13-4 

21 

0.17  ±  0.005 

49 

193 

0.10 

0.22 

13-5 

15 

0.18  ±  0.01 

32 

204 

0.09 

0.28 

13-1 

11 

0.23  ±  0  02 

47 

234 

0  17 

0.41 

13-6 

7.5 

032  ±  0.03 

% 

298 

0.35 

0.66 

13-7 

5.5 

0.53  ±  0.05 

73 

248 

0.45 

0.93 

13-8 

2.5 

0.75  ±  0.15 

40 

126 

0.53 

1.10 
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FIG.  S.  Total  surface  scattering  area  per  unit  volume  (5/V)  as  a 
function  of  densification  (i.e.,  percent  theoretical  density)  for  silica 
and  alumina. 


Stable  removal  of  pore  channels  one  by  one  from  an 
interconnected  network.  During  this  intermediate  stage, 
the  average  pore  channel  radius  remains  constant.  Only 
in  the  final  stage,  in  contrast  to  the  silica,  do  the  isolated 
pores  exist  and  require  the  proximity  of  grain  boundaries 
for  their  ultimate  removal. 


so  eo  70  eo  90  100 


%  TD 

FIG.  6.  Pore  volume  fraction,  normalized  to  5/V,  as  a  function  of 
densification  (%  TD)  for  silica  and  alumina. 


%  TD 


FIG.  7.  R,g(0)  as  a  function  of  percent  theoretical  density  for  silica 
and  alumina. 


VI.  SUMMARY  AND  CONCLUSIONS 

Small-angle  neutron  scattering  measurements  were 
used  to  examine  the  evolution  of  the  porous  micro- 
structure  of  ceramic  bodies  as  a  function  of  densification. 
Topologically,  both  systems  were  shown  to  evolve 
in  the  same  way.  However,  the  structural  evolution  in 
terms  of  pore  size  and  morphology  displays  a  different 
characteristic  signature  depending  on  the  sintering 
mechanism  that  is  operative.  For  the  densification 
of  a  glassy  ceramic,  porous  silica,  the  densification 
path  that  is  observed  is  attributed  to  the  viscous 
flow  mechanism  of  densification  and  to  the  absence 
of  a  grain  boundary.  For  the  densification  of  a  crystal¬ 
line  ceramic,  poly  crystal-line  alpha-alumina,  the  densifi¬ 
cation  path  that  is  observed  is  attributed  to  the 
dominance  of  difiiision  sintering  mechanisms  and 
to  the  conditions  imposed  by  the  existence  of  the 
grain  boundary.  Whereas  the  density  of  the  green 
body  established  the  radius  of  the  interconnected 
pore  channels,  which  was  constant  throughout  inter¬ 
mediate  stage  sintering  in  alumina,  the  density  of  the 
green  body  may  play  a  less  important  role  in  the  sinter¬ 
ing  of  glassy  silica. 
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